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Abstract 

This  manuscript  examines  transport  through  oxygen  ion  and  electronic  conducting  membranes  including  electrode/membrane  interfaces  by 
explicitly  incorporating  both  ionic  and  electronic  transport  through  the  membranes  and  across  interfaces.  Spatial  variation  of  electrochemical 
potential  of  oxygen  ions,  fc0 2- ,  electrochemical  potential  of  electrons,  /xe  (or  reduced  negative  electrochemical  potential  of  electrons,  <p  — 
—fie/e,  where  e  is  the  electronic  charge),  and  chemical  potential  of  molecular  oxygen,  /zo2,  through  membrane  and  across  interfaces  are 
examined  as  functions  of  transport  properties  of  membranes  and  interfacial  regions.  The  analysis  shows  that  description  of  transport  across 
electrode/membrane  interfaces  requires  two  transport  parameters — one  for  ions,  and  the  other  for  electrons.  The  transport  equations  are  applied 
to  fuel  cells,  pressure-driven  oxygen  separation  through  mixed  ionic-electronic  conducting  (MIEC)  membranes,  and  voltage-driven  oxygen 
separation  through  predominantly  oxygen  ion  conducting  membranes.  In  fuel  cells  and  MIEC  oxygen  separation  membranes,  the  /zo2  varies 
monotonically  between  the  two  end  values  corresponding  to  those  at  the  two  electrodes.  Thus,  in  fuel  cells  and  MIEC  oxygen  separation 
membranes,  the  stability  of  the  membrane  is  assured  as  long  as  the  oxygen  partial  pressure,  po2,  on  the  fuel  side  or  the  permeate  side  is 
above  the  decomposition  oxygen  partial  pressure  of  the  membrane.  By  contrast,  in  voltage-driven  oxygen  separation  membranes,  /xo2  in  the 
membrane  can  lie  outside  of  the  end  values.  Thus,  in  the  case  of  oxygen  separation  under  an  applied  voltage,  the  transport  properties  of  the 
material  and  the  interfaces  determine  membrane  stability.  Implications  of  the  analysis  concerning  the  applicability  of  the  so-called  the  three- 
electrode  system  under  an  applied  voltage  to  investigate  electrode  polarization  are  presented.  It  is  shown  that  the  use  of  the  three-electrode 
system  for  the  estimation  of  electrode  kinetics  can  lead  to  significant  errors  at  high  applied  voltages,  and  may  result  in  overestimation  of 
electrocatalytic  activity  of  the  electrode.  This  manuscript  also  defines  electrode  overpotential  in  terms  of  the  rate  of  potential  useful  work 
degraded  as  an  irreversible  process  at  an  electrode  and  the  net  current  measured  in  the  external  circuit.  The  as-defined  overpotential  may  not 
be  experimentally  measurable. 

©  2005  Elsevier  B.Y.  All  rights  reserved. 

Keywords:  Mixed  ion-electronic  conductors;  Oxygen  ion  conductors;  Fuel  cells;  Oxygen  separation;  Transport  across  interfaces 


1.  Introduction 

Many  oxides  transport  both  ionic  (cationic  and/or  anionic) 
and  electronic  species.  Of  particular  interest  for  applications 
such  as  fuel  cells,  water  electrolysis,  potentiometric  oxy¬ 
gen  sensors,  oxygen  separation  under  an  applied  voltage, 
etc.  are  materials,  which  can  predominantly  transport  oxygen 
ions.  For  oxygen  separation  from  air  under  an  applied  pres¬ 
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sure  gradient,  materials,  which  exhibit  significant  transport 
of  both  oxygen  ions  and  electronic  defects,  are  of  interest. 
Such  materials  are  known  as  mixed  ionic-electronic  conduc¬ 
tors  (MIEC).  Transport  in  oxygen  ion  conductors  or  in  MIEC 
materials  can  be  described  by  phenomenological  equations, 
which  have  been  well  established  for  over  50  years  [1-5]. 
Several  published  works  have  applied  these  transport  equa¬ 
tions  for  investigating  ionic  and  electronic  transport  through 
membranes  [6-11].  In  most  of  these  studies,  the  focus  was  on 
transport  through  membranes  with  emphasis  on  the  role  of 
defect  chemistry,  but  rarely  on  the  role  of  interfaces.  In  those 
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Nomenclature 

Bi 

mobility  of  species  i 

Ci 

concentration  of  species  i 

e 

electronic  charge 

E,E0 

Nernst  voltage 

Ea 

applied  voltage 

Em 

measured  voltage 

F 

Faraday  constant 

§e 

electronic  specific  conductance 

8i 

ionic  specific  conductance 

ji 

flux  of  species  i 

h 

electronic  current  density 

h 

current  density  due  to  species  i  (also  ionic  cur¬ 
rent  density) 

II 

load  current  (density) 

Boltzmann  constant 

i 

membrane  thickness 

«o2 

number  of  moles  of  O2  gas 

Po2 

partial  pressure  of  molecular  oxygen 

re 

electronic  area  specific  resistance 

n 

ionic  area  specific  resistance 

R 

ideal  gas  constant 

Re 

electronic  area  specific  resistance  of  the  mem¬ 
brane  including  interfaces 

Ri 

ionic  area  specific  resistance  of  the  membrane 
including  interfaces 

Rl 

load 

t 

time 

u 

ionic  transference  number 

T 

temperature 

f^pore 

volume  of  pore 

Zi 

valence  of  species  i 

Greek  letters 

5 

interface  thickness 

0 

overpotential 

Pi 

chemical  potential  of  species  i 

Pe 

electrochemical  potential  of  electrons 

Pi 

electrochemical  potential  of  species  i 

Po2 

chemical  potential  of  oxygen  gas 

Pq2- 

electrochemical  potential  of  oxygen  ions 

decomp 

Ro2 

decomposition  oxygen  chemical  potential  of 
the  membrane 

^membrane  chemjca]_  potential  of  oxygen  in  the  mem¬ 
brane 

A 

chemical  potential  of  oxygen  in  the  gas  phase 
close  to  the  electrode  I/membrane  interface 

Po2 

chemical  potential  of  oxygen  in  the  gas  phase 
close  to  the  electrode  II/membrane  interface 

Pe 

electronic  resistivity 

Pi 

ionic  resistivity 

°e 

electronic  conductivity 

Oi  conductivity  due  to  species  i  (also  ionic  con¬ 
ductivity) 

<p  reduced  (negative)  electrochemical  potential 

of  electrons 

cp1  reduced  (negative)  electrochemical  potential 

of  electrons  outside  the  membrane  close  to  the 
electrode  I/membrane  interface 
<p11  reduced  (negative)  electrochemical  potential 

of  electrons  outside  the  membrane  close  to  the 
electrode  II/membrane  interface 
0  electrostatic  potential  (Galvani) 


studies  which  have  examined  transport  through  membranes 
taking  into  account  interfaces,  the  focus  has  been  on  either 
ion  transport  across  interfaces  or  electron  transport  across 
interfaces,  but  usually  not  both.  Phenomenological  transport 
equations  for  the  flux  of  a  charged  species  are  usually  writ¬ 
ten  in  terms  of  the  gradient  in  electrochemical  potential  of 
the  species  involved  and  the  relevant  transport  parameters.  A 
fundamental  assumption  made,  albeit  often  tacitly,  is  that  of 
the  existence  of  local  equilibrium,  which  leads  to  relations 
between  electrochemical  potentials  of  charged  species  and 
chemical  potentials  of  neutral  species  [2,3].  For  example,  for 
the  case  of  oxygen  ions  and  neutral  oxygen  molecules  (or 
atoms),  local  equilibrium  may  be  described  in  terms  of  oxy¬ 
gen  molecules  (or  atoms),  electrons  and  oxygen  ions  as  the 
species  participating  in  the  reaction. 

The  above  general  framework  is  assumed  regardless  of 
the  prevailing  point  defects,  vacancies  or  interstitials  or  both. 
Additional  equations  can  be  written  in  terms  of  the  concentra¬ 
tions  and  transport  properties  of  point  defects.  This,  however, 
often  requires  one  to  make  simplifying  assumptions  concern¬ 
ing  the  thermodynamics  of  defects.  The  assumption  of  local 
equilibrium  is  universally  made,  either  explicitly  or  implic¬ 
itly  [12].  The  assumption  of  local  equilibrium  in  a  system 
also  implies  the  existence  of  equilibrium  at  any  position  in 
the  system,  wherein  the  system  may  be  single  phase  or  may 
consist  of  multiple  phases — including  the  gas  phase.  Trans¬ 
port  of  a  species  occurs  from  one  position  to  another  (ad¬ 
jacent)  position  in  response  to  difference  (gradient)  in  the 
relevant  thermodynamic  potential.  In  MIEC  materials,  the 
transport  of  both  ions  and  electronic  defects  is  taken  into  ac¬ 
count,  as  the  conductivities  of  the  two  species  are  nonzero 
finite,  and  often  comparable.  However,  when  dealing  with 
materials  that  are  predominantly  ionic  conductors  with  neg¬ 
ligible  electronic  conductivity,  often  it  has  been  the  prac¬ 
tice  to  ignore  the  electronic  transport  altogether  [13].  This 
may  be  satisfactory  if  the  interest  is  only  in  the  net  current 
flowing  through  the  material.  The  neglect  of  electronic  cur¬ 
rent,  however,  amounts  to  an  inconsistency  insofar  as  the 
assumption  of  local  equilibrium  is  concerned.  That  is,  the  as¬ 
sumption  of  local  equilibrium  and  the  simultaneous  neglect 
of  electronic  conduction  are  contradictory  assumptions.  If  lo- 
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cal  equilibrium  is  assumed,  which  is  almost  always  necessary 
for  writing  down  relevant  transport  equations,  the  electronic 
conductivity  cannot  be  assumed  to  be  identically  equal  to 
zero.  In  many  cases,  this  inconsistency  often  does  not  lead 
to  significant  problems.  However,  as  will  be  discussed  in  this 
manuscript,  and  has  been  discussed  to  some  extent  in  earlier 
manuscripts,  there  are  many  situations  in  which  the  neglect 
of  electronic  current,  however  small,  leads  to  substantial  dif¬ 
ficulties  in  analysis — and  thus  also  in  the  interpretation  of 
experimental  results  [14-16]. 

In  devices  such  as  fuel  cells  and  MIEC  oxygen  separa¬ 
tion  membranes,  no  external  voltage  is  applied.  Whatever 
voltage  is  developed  across  the  membrane  is  the  result  of 
differences  in  chemical  potentials  of  neutral  species  imposed 
and  the  transport  properties  of  the  membrane.  However,  when 
dealing  with  devices  such  as  oxygen  separation  membranes 
under  an  externally  applied  voltage,  the  system  then  has  an 
additional  (experimental)  degree  of  freedom,  which  has  im¬ 
portant  consequences  concerning  the  design  of  experiments 
and  the  interpretation  of  data,  especially  concerning  studies 
on  electrode  kinetics  or  in  oxygen  separation  under  an  applied 
voltage  or  water  electrolysis.  In  devices  such  as  solid  oxide 
fuel  cells,  it  is  generally  not  possible  to  obtain  an  accurate 
measurement  of  electrode  overpotential,  since  it  is  necessary 
to  place  a  reference  electrode  on  the  surface — unlike  in  aque¬ 
ous  electrochemistry  wherein  it  can  be  immersed  in  the  liquid 
electrolyte.  In  order  to  circumvent  this  problem,  often  the  so- 
called  three-electrode  system  is  used,  wherein  the  electrolyte 
is  in  the  form  of  a  thick  pellet.  The  working  electrode,  the 
counter  electrode,  and  the  reference  electrode  are  placed  in 
an  axisymmetric  arrangement.  Then,  to  investigate  electrode 
kinetics  at  the  working  electrode,  which  may  be  a  prospective 
cathode  for  a  fuel  cell,  an  external  voltage  is  applied  across  the 
counter  and  the  working  electrodes.  An  assumption  is  made 
that  the  externally  applied  voltage  is  essentially  equivalent 
to  that  generated  by  the  application  of  an  external  chemical 
potential  difference,  and  all  that  is  required  is  to  measure 
the  current  (density)  as  a  function  of  overpotential,  the  lat¬ 
ter  measured  with  respect  to  the  reference  electrode — and 
corrected  for  ohmic  losses.  This  approach  is  assumed  to  rep¬ 
resent  the  situation  at  a  real  cathode  in  a  fuel  cell.  As  will 
be  demonstrated  in  this  manuscript,  the  three-electrode  sys¬ 
tem  is  fundamentally  and  physically  different  from  a  fuel  cell 
and  its  use  under  an  externally  applied  voltage  to  investigate 
electrode  kinetics  can  lead  to  substantial  error — where  the 
magnitude  of  the  error  depends  upon  a  number  of  parame¬ 
ters,  including  the  magnitude  of  the  applied  voltage. 

This  manuscript  first  develops  transport  equations  us¬ 
ing  the  standard  phenomenological  approach.  The  difference 
compared  to  much  of  the  prior  work,  however,  is  the  in¬ 
clusion  of  the  electronic  current  in  a  predominantly  ionic 
conductor — and  especially  its  role  in  the  estimation  of  the 
local  chemical  potential  of  oxygen,  /x o2 ,  which  determines 
the  stability  of  the  material,  and  also  the  net  current  flow¬ 
ing  through  the  system.  This  also  implies  the  inclusion  of 
both  ionic  and  electronic  currents  through  the  membrane 


and  across  the  electrode/membrane  interfaces.  Three  cases  of 
practical  interest  are  examined:  (1)  oxygen  separation  using 
an  MIEC  membrane,  (2)  a  fuel  cell,  and  (3)  oxygen  separation 
under  the  application  of  an  external  voltage.  The  implications 
of  the  latter  are  examined  for  the  often-used  three-electrode 
system  under  an  applied  voltage  for  the  investigation  of  elec¬ 
trode  kinetics.  All  equations  are  developed  for  the  case  where 
the  transport  properties  of  a  given  region  are  constant,  inde¬ 
pendent  of  position  and  /xo2 .  Although  in  many  materials  the 
transport  properties  are  /xo2 -dependent — and  thus  are  func¬ 
tions  of  position — such  a  simplifying  assumption  is  deemed 
necessary  to  allow  for  the  development  of  simple  analyti¬ 
cal  equations.  This  approach  also  allows  for  description  in 
terms  of  simple  equivalent  circuits.  If  quantitative  relations 
between  transport  parameters  and  /x o2  are  known,  which  is 
rarely  the  case  for  most  materials  (and  interfacial  regions),  it 
would  be  a  straightforward  matter  to  numerically  solve  the 
relevant  transport  equations. 


2.  Theoretical  analysis 


In  what  follows,  transport  only  under  isothermal  con¬ 
ditions  is  addressed.  It  is  also  assumed  that  gas  transport 
through  the  porous  electrodes  to  the  electrode/membrane 
interfaces  is  sufficiently  fast  such  that  concentration  polar¬ 
ization  is  negligible.1  General  transport  equation  for  a  one¬ 
dimensional  flux  is  given  by  [2] : 


ji  =  -CiBiVfri  =  -CiBi 


dftj 

dx 


(assuming  the  flux  of  a  charged  species  occurs  down  its  elec¬ 
trochemical  potential  gradient)  in  #  (cm-2  s_1 ),  where  C/  and 
Bi  are,  respectively,  concentration  and  mobility  of  species  /, 
and  the  electrochemical  potential  of  species  i  is  given  by 


A/  =  +  Zie& 


where  /x/  is  the  chemical  potential  of  species  i  (J  or  ergs),  zt 
the  valence  of  species  i,  e  the  electronic  charge  (C),  and  0 
the  local  electrostatic  potential  or  Galvani  potential  (V). 

The  current  density  due  to  species  i  for  a  one-dimensional 
case  is  given  by 


<Ji 

h  =  zteji  = - VAi  = 

zie 


at  dfii 
Zte  dx 


in  A  cm-2,  where  0[  =  zfe2CjBi  is  the  conductivity  due  to 
species  i  (S  cm_1)  =  1/p;,  where  pi  is  the  resistivity  (Q  cm). 
In  the  above,  /x/  and  A/  are  defined  on  a  per  species  (per  ion, 
per  electron,  per  atom,  or  per  molecule)  basis.  If  defined  on 


1  In  practice,  this  assumption  often  may  not  be  valid.  For  example,  this 

assumption  may  not  be  valid  when  using  the  three-electrode  system  under 
an  externally  applied  voltage,  the  implications  of  which  are  discussed  later. 
Also,  this  assumption  is  almost  certainly  not  valid  in  devices  such  as  solid 
oxide  fuel  cells,  which  use  porous  electrodes  several  microns  in  thickness. 
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a  per  mole  basis,  e  is  replaced  by  the  Faraday  constant,  F 
(Cmor1). 

Local  equilibrium  is  assumed  to  prevail  everywhere  in  the 
system  [12].  The  cation  sublattice  is  assumed  to  be  rigid  in 
what  follows.  For  an  oxygen  ion  conductor,  the  equilibrium 
of  interest  is  the  following,  namely: 

i02  +  2e'  O2-  (4) 

Thus,  at  (local — at  any  position  in  the  system)  equilibrium 

[2,12] 

^/xo2  +  =  Iiqi-  (5) 

where 

A02-  =  fi02-  —  2  e<P  (6a) 

and 

fie  =  i±e  -  e@  (6b) 


Local  equilibrium  also  implies  that  any  changes  made  in 
chemical  or  electrochemical  potentials  (at  any  given  posi¬ 
tion)  must  satisfy  the  following  relation: 

o2  +  2  hjle  =  8/x02-  (7) 


where  8X  denotes  a  small  deviation  in  the  thermodynamic 
potential  X  from  the  equilibrium  state.  The  above  implies  that 
the  material  in  question  must  be  able  to  transport  at  least  two 
of  the  species  (at  least  two  of  O2  (O),  O2-  and  e ),  in  and  out, 
in  order  to  maintain  local  equilibrium.  This  also  implies  that 
when  potential  gradients  exist,  one  can  never  have  the  flux  of 
any  of  the  species  (participating  in  the  equilibrium  reaction) 
in  the  above,  identically  equal  to  zero.  That  is,  even  for  a 
predominantly  ionic  conductor,  we  must  always  have  Ie  /  0, 
even  though  it  may  be  small.  This  has  important  implications 
concerning  the  chemical  potential  of  oxygen  and  thus  the 
very  stability  of  the  material,  which  will  be  discussed  in  this 
manuscript. 

Ionic  and  electronic  current  densities  are  given  respec¬ 
tively  by2 


ao2~ 

2e 


Cfi  d  fa 
2e  dx 


where  i  is  written  for  O2  ,  and 


<7e 

V  fie  = 


e 


(Je  d fie 
e  dx 


Let  us  write 


fre 

(p  = - 

e 


fie 

e 


+  <P 


(9) 

(10) 


-02  +  2e^O~2 
2 

1 

~M02  +  2Pe  =M0-2 


Gas  Phase 


/ 


(p 


Membrane 


CX 


Pore 


-  02  +2e^O~2 
2 

1 

~Mo2  +  2Pe  =  P  0-2 


II 


Mo, 


(p 
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-On+2e^O~2 

2 

1 

-Mo2  +  2Me  =  M0-2 


Gas  Phase 


◄ -  l  - ► 


Fig.  1.  A  schematic  illustration  of  a  membrane  of  thickness  l,  which  can 
transport  O2-,  e,  and  thus  O,  across  which  exists  a  difference  in  chemi¬ 
cal  potentials  of  oxygen,  (p20  —  P1q0),  and  a  difference  in  reduced  (neg¬ 
ative)  electrochemical  potentials  of  electrons  ((pl  —  (p n).  Local  equilibrium 
in  the  system  implies  that  the  reaction  1 02  +  2er  — ►  O2-  is  in  equilib¬ 
rium  everywhere  in  the  system.  This  means  if  a  pore  of  volume  Vpore  ex¬ 
ists  at  some  position  in  the  membrane,  the  oxygen  partial  pressure  in  the 
pore  will  be  that  corresponding  to  the  local  chemical  potential  of  oxygen 
PO2  —  +  /?Tln  poo,  where  /Xq2  is  the  standard  state  chemical  poten¬ 

tial  and  the  ideal  gas  law  is  assumed.  Thus,  the  number  of  moles  of  02  in 
the  pore  is  given  by  no,  =  PCh  ^po re/RT.  If  now  Vp0re  — >  0,  then  no2  — >  0. 
Yet,  the  local  po2  (and  )  is  well  defined. 


which  is  the  Hebb  notation  [1].  Here  cp  is  the  (negative)  re¬ 
duced  electrochemical  potential  of  electrons  (V).  Using  this 
notation: 


dip 

dx 


(11) 


and 


Oi  dpio2  dip 

- &i - 

\e  dx  dx 


(12) 


Fig.  1  shows  a  schematic  illustration  of  a  membrane,  which 
can  transport  02~,e,  and  thus  O  (and  thus  also  O2),  at  the  two 
electrodes  of  which,  different  chemical  potentials  of  oxygen 
and  different  electrochemical  potentials  of  electrons  may  ex¬ 
ist.  It  is  assumed  that  local  equilibrium  prevails  at  all  positions 
in  the  system,  including  the  gas  phase  just  outside  the  mem¬ 
brane.  Fig.  1  also  shows  the  significance  of  local  equilibrium 
in  a  fully  dense  material,  via  an  example  of  a  pore  whose  vol¬ 
ume  approaches  zero.  If  a  small  pore  exists  at  any  position 
within  the  membrane,  the  local  /xo2  is  related  to  the  local 
oxygen  partial  pressure,  po2,  by  /xo2  =  /Xq?  +  ^Tln po2 , 
with  the  number  of  moles  of  O2  in  the  pore,  no2 ,  related  to 
po2 » the  pore  volume  and  temperature  by  the  ideal  gas  law.3 
If  the  pore  volume  approaches  zero,  the  number  of  moles  of 
oxygen  in  the  pore,  no2 ,  approaches  zero.  However,  the  /xo2 
and  the  po2  continue  to  exist. Transport  of  ions  and  electrons 
also  occurs  across  the  two-electrode/electrolyte  interfaces. 
The  interface  thicknesses  are  generally  not  known — and  not 
well  defined  either.  This  is  because  interfacial  regions  can 
often  be  diffuse — due  to  compositional  gradients  and  space 
charge  effects.  It  will,  nevertheless,  be  assumed  here  that  the 


2  Even  though  in  most  oxygen  ion  conductors,  transport  of  oxygen  ions 
occurs  by  a  vacancy  mechanism,  all  transport  equations  are  given  here  in  a 
generic  form — that  is,  without  reference  to  a  specific  defect  type. 


3  Or  another  appropriate  gas  law,  depending  upon  the  pressure  and  tem¬ 
perature.  In  this  manuscript,  it  is  assumed  that  the  ideal  gas  law  is  applicable 
over  the  entire  pon  range. 
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‘interfacial’  regions  are  of  certain  thicknesses  and  describe 
transport  properties  through  the  two  interfaces  in  terms  of  the 
transport  properties  of  interfacial  regions  of  certain  properties 
and  thicknesses: 

cr-  and  8 ' — ionic  conductivity  and  thickness  of  electrode 
I/membrane  interface; 

o'e  and  8' — electronic  conductivity  and  thickness  of  elec¬ 
trode  I/membrane  interface; 

cr "  and  8" — ionic  conductivity  and  thickness  of  electrode 
II/membrane  interface; 

cr"  and  8" — electronic  conductivity  and  thickness  of  elec¬ 
trode  II/membrane  interface. 


possible  to  separately  determine  conductivities  and  thick¬ 
nesses  of  interface  regions,  the  specific  interface  conduc¬ 
tances  (gj,  g",  g'e  andg")  or  the  specific  interface  resistances 
(r-,  r",  r^andr")  can  in  principle  be  measured  experimen¬ 
tally. 

Using  the  above  terminology,  ionic  and  electronic  current 
densities  across  the  interfaces  are  given  by5 

7  =  -|4/4>2  -  m o2)  +  SiW1  -  <p') 

=  +  (i5) 

4e  ri  ri 


Thus,  specific  ionic  conductance  of  interface  I,  namely  g-,  is 


(13a) 


in  £2_1  cm-2  or  S  cm-2,  and  r\  is  the  specific  ionic  resistance 
of  interface  I  in  ^  cm2,  and  specific  electronic  conductance 
of  interface  I,  namely  gre,  is 


(13b) 


in  £2_1  cm-2  or  S  cm-2,  and  rJe  is  the  specific  electronic  re¬ 
sistance  of  interface  I  in  Q  cm2.  Similarly,  specific  ionic  and 
electronic  conductances  of  interface  II  are  given  by 


U! 

G_j_ 

8" 


and 


a 


a 

e 


8” 


(13c) 


(13d) 


respectively.  The  preceding  assumes,  for  simplicity,  that  the 
transport  properties  over  the  interface  thickness  are  constant. 
This  of  course  need  not  be  the  case,  and  general  equations 
can  be  written  as  follows: 


(14) 


and  similarly  for  the  other  specific  conductances  or  resis¬ 
tances.  The  spatial  dependence  of  transport  properties  im¬ 
plied  in  Eq.  (14)  is  through  their  dependence  on  po2.  In 
what  follows,  transport  properties  of  the  interfacial  regions 
are  assumed  to  be  constant  over  the  thickness  and  indepen¬ 
dent  of  po2 .4  It  will  also  be  assumed  that  the  interface  thick¬ 
nesses  are  much  smaller  than  the  membrane  thickness,  that 
is,  8\  8"  <£.  It  is  to  be  noted  that  while  it  is  generally  not 


4  In  general,  the  transport  parameters  of  interfacial  regions  are  functions  of 
po2 .  For  example,  a  simple  analysis  based  on  adsorption  shows  that  usually 
r ■  oc  pQlJn  ,  where  n'  is  a  positive  number  typically  an  integer.  This  aspect  is 
ignored  here  since  the  objective  is  to  provide  analytical  equations  to  describe 
general  features.  For  any  given  material,  if  the  dependence  is  known,  it  is  a 
trivial  matter  to  incorporate  it  and  solve  equations  numerically. 


and 


a” 

6(  f  //  II  l  i  ft  t  ft  Hi 

~ r-o2 )  +  Si{<P  ~<p  ) 
1  Mo2  _  r02  <p  -  <pn 

~4e  V!  +^T~ 


(16) 


as  the  ionic  current  densities  across  the  two  interfaces,  and 


(17) 


(18) 


as  the  electronic  current  densities  across  the  two  interfaces.  In 
the  preceding  equations,  p'0i  and  <pf  are,  respectively,  chem¬ 
ical  potential  of  oxygen  and  reduced  (negative)  electrochem¬ 
ical  potential  of  electrons  in  the  membrane,  just  inside  elec¬ 
trode  I/membrane  interface;  /Xq2  and  <pn  the  corresponding 
quantities  in  the  membrane  just  inside  interface  II.  In  steady 
state:6 


and 


(19) 


(20) 


where  //  and  Ie  are,  respectively,  ionic  and  electronic  current 
densities  through  the  bulk  of  the  membrane.  In  general,  trans¬ 
port  properties  of  most  electrolyte  or  MIEC  materials  depend 
upon  the  local  oxygen  partial  pressure,  po2 ,  or  chemical  po¬ 
tential,  po2  •  Thus,  the  relevant  equations  are  (11)  and  (12), 
where  the  transport  properties  are  po2  -dependent,  and  thus 
position  dependent.  Integration  of  Eq.  (12)  gives  [3] 


1 

4e 


jt 


(7i  d/io: 


(21) 


5  If  the  transport  properties  vary  across  the  interfaces,  then  thq  relevant 
transport  equations  are  of  the  form  /;  =  2  |  ^  fQ  cr-  d/xcb  —  fQ  cr-  dip  j . 

6  In  three  dimensions,  the  criterion  for  steady  state  is  given  by  V  f  =  0  and 
V  • Ie  =  0,  where  7/  and  7e  are  vectors. 
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and  integration  of  Eq.  (11)  gives 


(22) 


where  the  integration  is  between  the  two  (just  inside)  elec¬ 
trode/membrane  interfaces.  In  what  follows,  ionic  and  elec¬ 
tronic  conductivities  of  the  membrane  will  be  assumed  to  be 
constant,  independent  of  position.  This  allows  explicit  deter¬ 
mination  of  fluxes  in  terms  of  transport  properties,  and  more 
importantly  the  estimation  of  /Xqo  and  /Xq2  ,  and  cp'  and  cp"  in 
terms  of  measurable  parameters.  Also,  this  makes  it  possi¬ 
ble  to  analyze  a  number  of  cases  of  interest  analytically.  The 
preceding  equations  then  become 


4  el 


{(Aq2  -  Mo2)  + 


(23) 


and 

Ie  =  j{cp'-(p"}  (24) 


The  parameters,  which  are  known  (or  which  can  in  principle 
be  determined  experimentally)  are  the  following:  <7;,  cre,  g-  = 
(!/?•,■),  g'e  =  (1/r'X  g'l  =  (1  /r"),  g'e  =  (1  /r'e),  C  Po2,  II o2, 
<pl,  cp n,  Ii,  and  Ie.  The  unknowns  are:  /Xq2,  /Xq2,  cp\  and  cp" . 

From  the  transport  of  electrons  across  interface  I: 


/  I  Le  I  /  T 

<P  =(P - 7  =(P  ~reIe 


g 


e 


From  the  transport  of  electrons  across  interface  II: 


(25) 


<p”  =  A1  + 


g 


£  _  II  I  //  j 

-=(p  +reIe 

e 


(26) 


From  the  transport  of  oxygen  ions  across  interface  I: 


=  Pq2  +  4e(r'li  -  r'eIe ) 

(27) 


From  the  transport  of  oxygen  ions  across  interface  II: 


//  II  A 

Mo2  =  ^o2  -  4e 


a”  l-  —  o'!  J 
c>e  c>i  1e 


g'eg'i 


Eqs.  (25)-(28)  describe  the  unknown  potentials,  namely,  /Xq2  , 
/Xq2,  <p’ ,  and  ip" ,  in  terms  of  transport  properties  of  interfaces 
and  the  net  ionic  and  electronic  current  densities,  namely  f 
and  Ie.  An  important  point  to  note  from  Eqs.  (27)  and  (28)  is 
that  even  for  a  predominantly  oxygen  ion  conducting  mem¬ 
brane,  one  cannot  a  priori  assume  either  r'e  or  r"  to  be  infinite, 
which  makes  Ie  zero  and  the  fio2  in  the  membrane  indetermi¬ 
nate  (product  of  oo  and  0)  and  violates  the  criterion  of  local 
equilibrium.7.  Thus,  r'e  and  r"  can  be  very  large,  but  not  in¬ 
finite.  Eqs.  (27)  and  (28)  also  show  that  usually  /Xq2  ^  /Xq2 

and  /Xq2  7^  /Xq2  ,  as  this  would  require  terms  in  parentheses 


7  Obviously,  re  cannot  be  infinite  either,  as  discussed  in  [14]. 


to  be  identically  zero — which  would  a  very  rare  and  a  special 
case,  and  not  a  general  case.  As  evident  from  Eqs.  (27)  and 
(28),  the  chemical  potentials  of  O2  just  inside  the  interfaces 
(in  the  membrane)  depend  upon  the  values  of  r-,  r'e ,  r",  r",  Ii , 
and  Ie. 

In  terms  of  the  reduced  electrochemical  potentials  of  elec¬ 
trons  at  the  two  electrodes,  the  net  electronic  current  density 
is  given  by 


Ie  = 


<p1-<p 


11 


lpe  +r'+r 


a 


(29) 


and  in  terms  of  the  chemical  potentials  of  oxygen  and  the 
reduced  electrochemical  potentials  of  electrons  at  the  two 
electrodes,  the  net  ionic  current  density  is  given  by 


1  Mq2  -  Pp2  (pl  -  A 

4e  Ipi  +  r ■  +  r"  Ipt  +  r-  +  r " 


(30) 


Eqs.  (29)  and  (30)  for  the  electronic  and  ionic  current  densi¬ 
ties,  respectively,  are  given  in  terms  of  parameters  which,  in 
principle,  are  measurable.  In  what  follows,  the  above  equa¬ 
tions  are  used  to  examine  the  various  cases  outlined  earlier. 


2.7.  The  generalized  form  of  chemical  potential  of 
oxygen  as  a  function  of  position  in  a  membrane  without 
the  incorporation  of  interfaces 


In  what  follows,  let  us  first  examine  the  chemical  po¬ 
tential  of  oxygen  as  a  function  of  position  in  a  membrane, 
without  the  incorporation  of  the  effect  of  interfaces.  This  ap¬ 
proach  leads  to  simple  analytical  expressions.  Let  us  consider 
a  membrane  of  thickness  f.  The  chemical  potentials  of  oxy¬ 
gen  at  the  two  sides  are  /Xq2  at  v  =  0  and  /Xq2  at  x  =  l.  The  net 
electronic  and  ionic  current  densities  through  the  membrane 
are  given  by  Eqs.  (11)  and  (12).  In  general,  the  conductiv¬ 
ities  are  functions  of  position  (which  actually  are  functions 
of  position-dependent  chemical  potential).  Summing  the  two 
current  densities,  rearranging  the  terms  and  integrating  gives: 


Mo2  (x)  =  p02  +  4e(li  +  Ie) 


dx 

Oifx) 


+  4e 


afx)  +  ae(x) 
&i(x) 


d  cp(x) 


(31) 


The  preceding  shows  that  the  chemical  potential  of  oxy¬ 
gen  within  the  membrane  can  be  completely  described  by 
Eq.  (31)  provided  the  net  current  flowing  through  the  mem¬ 
brane,  namely  7/  +  Ie  is  known,  and  transport  properties  are 
known  as  functions  of  position.  Eq.  (3 1)  is  applicable  for  both 
steady  and  transient  states.  In  the  transient  state,  all  depen¬ 
dent  variables  are  also  time  dependent,  that  is  crfx,  t ),  cre(x , 
t),  pio2(x ,  t ),  <p(x,  0,  7z(x,  t)  and  Ie(x,  t).  Yet,  Kirchoff’s  laws 
are  assumed  applicable,  which  gives  Ifx,  t)  +  Ie(x,  t)  =  7totai(0- 
That  is,  the  net  current  density,  7total(0>  is  a  function  of  time 
alone.  In  what  follows,  we  will  only  consider  the  steady  state. 
In  the  preceding  it  is  assumed  that  the  chemical  potentials  of 
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Fig.  2.  A  schematic  illustration  of  an  MIEC  membrane  under  an  applied  dif¬ 
ference  of  chemical  potential  of  oxygen,  with  no  applied  voltage  difference. 
Electronic  and  ionic  currents  are  of  equal  magnitude  but  are  in  the  opposite 
directions. 

oxygen  just  inside  the  two  interfaces  are  the  same  as  in  the  lo¬ 
cal  atmospheres.  The  exact  details  of  the  variation  of  fio2(x) 
with  position  depends  upon  the  magnitudes  and  spatial  de¬ 
pendencies  of  the  transport  properties  of  the  membrane,  the 
end  values  of  the  chemical  potentials,  namely  /Xqo  and  /Xq2, 
and  the  end  values  of  the  reduced  electrochemical  potentials 
of  electrons,  namely  cp1  and  cp11.  If  analytical  relations  are 
known  between  transport  properties  and  position,  Eq.  (31) 
can  be  readily  used  to  determine  fio2(x)  as  a  function  of 
position.8 * *  This  is  a  special  case  wherein  an  assumption  is 
made  that  /Xq2  =  and  /Xq2  =  /Xq0  .  However,  there  is  no 
justification  for  this  assumption  in  a  general  case,  as  already 
stated  previously.  That  is,  in  a  general  case,  /Xq2  ^  /Xq2  and 

/Xq2  7^  /Xqo  .  Under  such  conditions,  it  is  not  possible  to  ob¬ 
tain  a  simple  analytical  equation  (such  as  the  one  given  in  Eq. 
(31))  with  transport  properties  dependent  upon  position.  The 
objective  of  this  work  is  to  explore  the  role  of  interfaces  and 
obtain  explicit  expressions  for  /Xqo  and  /Xqo  .  Hence,  a  sim¬ 
plifying  assumption  is  made  here  that  transport  parameters 
are  spatially  invariant  in  the  membrane.  In  order  to  account 
for  the  interfaces,  as  stated  earlier,  properties  of  the  interfaces 
are  separately  included. 

2.2.  An  MIEC  oxygen  separation  membrane 

In  what  follows,  the  transport  properties  are  assumed  to 
be  independent  of  position  within  the  membrane.  In  this  case 
of  an  MIEC  oxygen  separation  membrane,  no  voltage  is  ap¬ 
plied  across  the  membrane.  However,  a  chemical  potential 
difference  of  oxygen  is  imposed  across  the  membrane.  As 
there  is  a  chemical  potential  difference  across  the  membrane, 
a  corresponding  difference  in  electrochemical  potential  of 
electrons  is  created.  Fig.  2  shows  a  schematic  illustration  of 
a  MIEC  membrane  across  which  there  is  a  gradient  in  the 
chemical  potential  of  oxygen.  Let  /Xqo  >  /Xq0  .  Since  neither 
an  external  voltage  is  applied,  nor  the  membrane  is  externally 
shorted,  the  net  current  is  zero.  That  is,  7Z  +  Ie  =  0  or  Ie  =  — /;. 

For  /Xq2  >  /Xq2  ,  oxygen  ion  transport  occurs  from  I  to  II 
(from  left  to  right  in  Fig.  2).  That  is  //  <  0,  and  therefore,  Ie  >  0. 


8  An  important  point  to  note  is  that  the  transport  properties  are  functions 

of  iio2  >  and  their  dependence  on  position  is  through  the  dependence  of  the 

chemical  potential  on  position. 


n 


<P 


—// 
Mo-* 


Fig.  3.  Schematic  variations  of  /zo2,  P02-  and  (p  for  the  case  depicted  in 
Fig.  2.  Abrupt  changes  in  /xo2,  (p,  and  jl02-  are  shown  across  interfaces 
since  8',  8"  <&  l . 


The  chemical  potential  of  oxygen  just  below  the  interface  at 
I  given  by  Eq.  (27)  becomes 


=  Mo2  -  Mr-  +  r'e)\Ii 


(32) 


That  is,  /Xq2  <  Mo2  • 

Similarly,  the  reduced  (negative)  electrochemical  poten¬ 
tial  of  electrons  just  inside  interface  I  is  given  by 


e 


8 


=  <P  + 


e 


g 


e 


(33) 


That  is,  cpf  <cpl. 

Similar  equations  can  be  written  for  /xo2  and  <p  just  inside 
interface  at  II,  namely  /Zq2  and  cp" .  For  example  note  that 


That  is,  Mq2  >  Pf  • 

The  reduced  (negative)  electrochemical  potential  of  elec¬ 
trons  just  inside  interface  II  is  similarly  given  by 


(35) 


That  is,  <p"  >  (p11. 

Note  also  that  Ie  =  —  <p''}  >  0  or  cp'  ><p" .  Rear¬ 

rangement  of  Eq.  (23)  gives 


4  ef 

{Mo2  -  Mo2i  = - h  +  W  -  <p"}  (36) 

Oi 


Since  //  <  0  and  <pf  >  (p" ,  it  is  clear  that  /Xq2  >  /Xq?  . 

Thus,  for  an  MIEC  membrane,  for  the  experimentally  se¬ 
lected  conditions  corresponding  to  /Xqo  >  /Xq2  ,  it  is  seen  that 

l4)2  >  n'o2  >  l4,2  >  llo2  and  <pl  >  <p'  >  <p"  >  ip11,  regard- 
less  of  the  absolute  values  of  transport  parameters  of  the  bulk 
membrane  and  interfaces.  Fig.  3  shows  the  schematic  varia¬ 
tions  of  iio2  and  <p  for  this  case,  where  it  is  tacitly  assumed  that 
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re  ~  re  ~  rc-  Finally,  as  the  transport  of  oxygen  ions  occurs 
from  left  to  right,  note  that  Aq2-  >  Aq2-  >  Aq2-  >  Aq2-  • 
The  spatial  variation  of  Aq2-  is  aiso  shown  in  Fig.  3. 

Since  p0i  >  n'G_,  >  >  pf ,  as  long  as  /iff  > 

li02  ,  where  /Xq2  is  the  chemical  potential  or  oxygen 
below  which  membrane  decomposition  can  occur,  the  stabil¬ 
ity  of  the  membrane  is  assured.  An  important  point  to  note, 
however,  is  that  the  membrane  may  continue  to  be  stable 

even  if  <  /Xq2  ,  provided  >  fi0i  .  That  is,  in¬ 
terface  properties  are  expected  to  have  a  profound  effect  in 
dictating  membrane  stability.  The  interface  properties  may 
also  depend  upon  the  atmosphere.  For  example,  for  a  given 
/Xq2  (and  thus  for  a  given  Pq2),  the  may  be  different 
in  H2/H2O  and  CO/CO2  atmospheres.  Thus,  it  is  possible 
that  a  given  MIEC  may  be  stable  (permeate  side)  in  one  gas 
mixture  but  not  in  the  other,  even  though  the  po2  in  the  two 
atmospheres  is  identical. 

Substitution  for  /Xq2  ,  /Xq2  ,  <pf,  and  ip"  in  terms  of  respective 
current  densities  and  transport  properties  of  the  two  interfaces 
into  Eq.  (23)  gives 


11 


i 

—  + 
<*i 


1 

Ae 


Mq2}  "F  +  ri  +  rfl 


(37) 


Also,  from  the  equations  for  7/  and  Ie ,  and  noting  that 
Ii  +  Ie  =  0,  using  Eqs.  (29)  and  (30)  it  is  easy  to  see  that 


T 


II 


ip  -ip  = 


ipe  +  A  +  r 


// 


A 


/A1 
^ o2 


[lpe  +  r'e  +  r"]  +  \ipi  +  r\  +  r"]  Ae 


(38) 


Thus,  in  general 


I  II  ^02  ^O: 

ip  -  ip  <  — - - 

Ae 


If  [lpe  +  re  +  re]f>  [ipi  +  r\  +  r"],  then 


CJ  II  ^  ^02  ^O: 

ip  -ip  ~ 


=  E0 


(39) 


where  Eq  is  the  Nernst  potential. 

Equivalent  circuit.  It  is  customary  to  represent  an  equiv¬ 
alent  circuit  for  a  bulk  MIEC  in  terms  of  a  Nernst  voltage, 
generated  by  the  imposed  chemical  potential  difference,  and 
ionic  and  electronic  resistances  [17].  However,  for  interfaces, 
the  common  practice  has  often  been  to  use  a  capacitor  in 
combination  with  a  resistor  [18].  The  charged  capacitor  then 
models  the  voltage  difference  created  by  the  chemical  poten¬ 
tial  difference,  and  in  fact  a  so-called  chemical  capacitance 
can  be  defined,  which  is  meant  to  embody  the  relevant  ther¬ 
modynamics.  The  concept  of  chemical  capacitance  has  also 
been  used  by  some  to  model  a  bulk  MIEC  membrane.  In  such 
a  case,  a  (chemical)  capacitor  is  introduced  between  two  sets 
of  parallel  ionic  and  electronic  resistor  segments.  Fig.  4(a) 
shows  a  schematic  illustration.  In  steady  state,  no  current 


flows  through  the  capacitor,  while  the  ionic  current  is  uniform 
in  the  ionic  segments,  and  the  electronic  current  is  uniform 
in  the  electronic  segments.  In  a  transient  state,  the  ionic  cur¬ 
rents  in  the  two  ionic  segments  are  different  and  so  are  the 
electronic  currents  in  the  two  electronic  segments.  However, 
Kirchoff ’s  laws  continue  to  hold,  that  is,  V  •(/;  +  Ie )  =  0,  which 
for  a  one-dimensional  case  reduces  to  7Z  +Ie  =  constant.  Also, 
in  the  transient  state,  there  is  a  net  current  flowing  through 
the  capacitor. 

In  this  manuscript,  however,  internal  EMF  sources  will 
be  used  to  model  bulk  as  well  as  interfacial  regions.  Since 
the  chemical  potential  varies  spatially  (including  across  inter¬ 
faces),  representation  of  interfacial  regions  in  terms  of  Nernst 
voltages  is  a  natural  extension  of  the  usual  practice  of  rep¬ 
resenting  bulk  regions  in  terms  of  Nernst  voltages.  That  is, 
there  is  no  particular  requirement,  which  states  that  one  must 
use  capacitors  to  describe  the  interfacial  regions,  if  an  EMF 
source  is  used  to  describe  the  bulk  regions. 

In  what  follows,  rationale  for  the  choice  of  internal  EMFs 
to  describe  both  bulk  and  interfacial  regions  of  an  MIEC 
under  an  applied  chemical  potential  is  first  described.  Con¬ 
sider  the  application  of  a  net  chemical  potential  differ¬ 
ence  across  a  membrane  such  that  the  position-dependent 
chemical  potential  of  oxygen  given  by  fio2(x ),  is  a  con¬ 
tinuous  function  of  position,  v.  Thus,  the  chemical  poten¬ 
tial  at  x  +  dx  is  given  by  pio2(x  +  dv)  =  /jlo2(x)  +  d/xo2(v). 
The  corresponding  Nernst  voltage  across  the  slice  dx  is 

given  by  d E{x)  =  dll°^Ax\  iocai  ionic  and  electronic 
area  specific  resistances  for  the  slice  are  dr;(v)  =  pi(x)  dx  and 
dre{x)  =  pe{x)  dx.  The  corresponding  equivalent  circuit  for  the 
slice  is  shown  in  Fig.  4(b).  Note  that  the  electronic  and  the 
ionic  segments  are  joined  at  the  ends  (nodes).  The  signifi¬ 
cance  is  that  the  existence  of  local  equilibrium  correspond¬ 
ing  to  the  reaction  5O2  +  2ef  02~  requires  the  transport 
of  ionic  and  electronic  currents  to  any  position,  which  corre¬ 
sponds  to  joining  of  the  electronic  and  the  ionic  segments  at 
the  nodes.  The  equivalent  circuit  for  the  entire  membrane  (by 
discretizing  it  in  slices  of  thickness  dx)  is  shown  in  Fig.  4(c). 
Now,  no  charge  accumulation  can  occur  at  any  node.  In  such  a 
case,  Kirchoff ’s  equations  must  be  applicable,  which  means 

7/(x,  t)  T  7^(v,  t)  =  Ii(x  T  dv,  t)  H-  7^(v  T  dv,  t)  —  •  *  • 

=  Ii(x  +  n  dx ,  t)  +  Ie{x  +  n  dx,  t)  (40) 

where  t  denotes  time  and  n  the  number  of  differential  seg¬ 
ments.  Eq.  (40)  is  applicable  both  under  transient  and  steady 
state  conditions.  For  this  reason,  the  time  dependence  is  ex¬ 
plicitly  included.  Under  transient  conditions,  ionic  current 
entering  a  node  may  be  different  from  that  leaving  the  node, 
and  similarly  for  the  electronic  current,  yet  consistent  with 
the  above  identity.9  In  such  a  case,  all  dependent  variables, 


9  The  same  situation  exists  if  one  uses  a  capacitor  (Fig.  4(a))  to  describe 
transport.  In  this  case,  the  magnitude  of  the  difference  in  ionic  currents  en¬ 
tering  and  exiting  the  capacitor  is  the  same  as  the  magnitude  of  the  electronic 
currents  entering  and  exiting  the  capacitor. 
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Fig.  4.  (a)  An  equivalent  circuit  describing  spatial  variation  of  transport  parameters  through  an  MIEC  using  a  chemical  capacitance,  (b)  An  equivalent  circuit 
describing  spatial  variation  of  transport  parameters  through  a  slice  dx  of  an  MIEC  using  an  internal  EMF  source,  which  is  related  to  local  variation  in  chemical 
potential  of  oxygen,  d/xo,  Or),  (c)  An  equivalent  circuit  describing  spatial  variation  of  transport  parameters  through  an  MIEC.  The  equivalent  circuit  is  applicable 
for  both  transient  and  steady  states.  For  the  transient  state,  the  ionic  and  electronic  currents  are  functions  of  position  and  time.  However,  Kirchoff ’s  laws  continue 
to  be  valid.  In  steady  state,  the  ionic  and  electronic  currents  are  independent  of  time  and  position  (for  a  one-dimensional  problem),  (d)  An  equivalent  circuit 
describing  spatial  variation  of  transport  parameters  through  an  MIEC.  The  equivalent  circuit  is  applicable  only  for  steady  state. 


including  /zo2C*),  are  time  dependent.  Thus,  in  a  transient 
state,  Ii(x,  t)  +  Ie(x ,  t)  =  Ii(x  +  d x,t)  +  Ie(x  +  dr,  t),  however //(x, 
t)  zfz  Ii(x  +  dx,  t)  and  4(x,  t)  ^  I e(x  +  dx,  t).  The  corresponding 
equivalent  circuit  is  shown  in  Fig.  4(c).  In  such  a  case,  the 
electronic  and  the  ionic  segments  must  be  joined  at  the  nodes. 
Also,  in  such  a  case,  the  ionic  and  the  electronic  current  den¬ 
sities  are  in  general  time  dependent. 

In  steady  state,  in  addition,  the  divergence  of  ionic 
and  electronic  currents  must  be  individually  zero,  that  is, 
V  •//  =  0  and  V-U  =  0,  and  obviously  current  densities  are 
not  time  dependent.  For  a  one-dimensional  case,  this  im¬ 
plies  that  Ii(x)  =  Ii(x  +  dx)  =  Ii(x  +  2  dx)  =  •  •  •  =  /z-(x  + 
n  dx)  =  •  •  •  and  Ie{x)  =  Ie(x  +  dx)  =  Ie(x  +  2  dx)  =  •  •  •  = 
Ie(x  +  n  dx)  =  •  •  -.Under  such  conditions,  the  equivalent  cir¬ 
cuit  may  also  be  described  by  that  in  Fig.  4(d).  Neverthe¬ 
less,  it  is  to  be  recognized  that  the  general  equivalent  cir¬ 
cuit  applicable  for  both  the  transient  and  steady  states  is 
that  given  in  Fig.  4(c).  If  this  equivalent  circuit  is  used  for 
steady  state,  then  it  should  be  further  specified  that  both  the 
ionic  and  electronic  currents  satisfy  the  conditions  V  •//  =  0 
and  V-/e  =  0.  For  a  one-dimensional  case,  this  is  equiva¬ 
lent  to  ionic  and  electronic  currents  being  spatially  invariant. 


If  one  uses  a  capacitor  instead,  then  the  equivalent  circuit 
given  in  Fig.  4(a)  is  applicable  in  both  transient  and  steady 
states.  For  the  steady  state,  no  current  flows  through  the 
capacitor. 

The  choice  of  /xo2 -independent  and  position-independent 
properties  for  any  given  region  allows  for  a  simple  equivalent 
circuit  representation  of  the  MIEC  membrane  just  described. 
The  corresponding  equivalent  circuit  is  given  in  Fig.  5.  Note 
again  the  representation  in  terms  of  internal  EMF  sources 
rather  than  charged  capacitors — which  is  consistent  with 
the  spatial  variation  in  /xo2.  The  internal  EMFs  are  given 
by 


(41) 


E,  = 


w 


4e 


RT 
4 F 


In 


W, 

- In 

4e 


(42) 
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MIEC  Membrane 


Fig.  5.  An  equivalent  circuit  for  the  MIEC  membrane.  The  circuit  elements  between  the  two  filled  elliptical  symbols  are  not  physically  separable;  they 
collectively  represent  the  equivalent  circuit.  All  E' s  (E' ,  E")  are  of  the  same  sign. 
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membrane  material,  the  latter  given  by 
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'  n  +  re 


(50) 


and  their  sum  is  given  by 


E=E'  +  Ei  +  E  = 


n  A  ~  Ao2 


\e 


V  A, 


RT 
4 ~F 


In 


(44) 


The  equivalent  circuit  can  be  readily  used  to  solve  for  the 
currents  flowing  through  the  membrane,  subject  to  the  con¬ 
dition  that  in  steady  state  /■  =  /"  =  f  and  /'  =  I'J  =  Ie .10 
And,  also  in  this  case  of  an  MIEC  membrane,  f  +  Ie  =  0.  Let 
us  write 


Ri  =  r\  +  tpi  +  r"  =  r\  +  rt  +  r'(  (45) 

and 

Re  =  r'e+  ipe  +  r"  =  re  +  re  +  r ''  (46) 

where  /y  =  lpt  =  Hot  and  re=lpe  =  Uoe.  It  is  easy  to  see  that 


E 

Ri  +  Re 


(47) 


The  measured  voltage  across  the  MIEC  membrane  is 


EM  =  <fl  -  P  =  leRe  = 


ER 


e 


Ri  +  Re 


(48) 


Thus,  the  effective  ionic  transference  number  of  the  mem¬ 
brane,  including  interfacial  effects,  is  given  by 


Em  Re  _ r’e  +  re  +  r" _ 

E  ~  Ri  +  Re~  r[  +  n  +  r’!  +  r'  +  re  +  r'J 


(49) 


Often,  the  measurement  of  open  circuit  voltage  across  an 
MIEC  is  used  to  estimate  the  transference  number  of  the 


10  The  equivalent  circuit,  as  stated  above,  can  also  be  used  for  analysis  of 
the  transient  problem  wherein  V  f  ■=/=■  0  and  V  • le  0,  although  V  (/z  +  Ie)  =  0, 
as  discussed  for  one-dimensional  cases  in  [15,16]. 


It  is  clear  that  only  in  a  special  case  when  re  y>>  rfe  +  r"  and 
ri  +  re  r-  +  r"  +  r’e  +  rne,  will  the  measurement  of  open 
circuit  voltage  (OCV)  give  the  material  ionic  transference 
number,  tf1.  That  is,  in  general,  the  measurement  of  OCV 
cannot  be  used  to  determine  the  ionic  or  the  electronic  trans¬ 
ference  numbers  of  the  material.  This  conclusion  is  similar  to 
the  one  arrived  at  by  Liu  and  Hu,  who  evaluated  the  effect  of 
interfacial  polarization  on  the  apparent  transference  numbers 
of  MIEC  materials  [19].  The  effect  of  interfacial  resistance 
on  the  apparent  transference  numbers  in  MIEC  materials  has 
also  been  examined  by  Kharton  and  Marques  [20] . 

2.3.  A  solid  oxide  fuel  cell 

In  a  solid  oxide  fuel  cell,  one  electrode  is  exposed  to  fuel 
(anode)  and  the  other  electrode  is  exposed  to  oxidant  (cath¬ 
ode).  The  membrane  is  of  a  predominantly  ion  conducting 
material,  usually  an  oxygen  ion  conductor.  During  opera¬ 
tion,  an  external  load  is  connected  such  that  a  net,  finite, 
nonzero  current  flows  through  the  membrane  (and  of  course 
through  the  external  circuit).  Thus,  the  net  current  is  not  zero, 
that  is  f  +  Ie^  0.  In  what  follows,  electrode  I  is  the  cathode, 
and  electrode  II  is  the  anode.  During  cell  operation,  oxy¬ 
gen  ions  transport  from  I  to  II  through  the  membrane.  It 
is  also  clear  that  /Xqo  >  /Xq0.  Note  also  that  (pl>cp n.  Since 
the  membrane  is  predominantly  an  ionic  conductor,  that  is, 
re  +  r'e  +  r"e  >  rt  +  r\  +  r",  it  is  clear  that  |/;|  >  \Ie\-  It  is  to 
be  emphasized  that,  for  the  reasons  already  discussed,  there 
must  be  a  finite,  nonzero  Ie ,  to  ensure  that  local  equilibrium 
is  established.  Fig.  6  shows  a  schematic  illustration. 

In  this  case  also,  Eq.  (37),  which  gives  the  difference  be¬ 
tween  the  reduced  electrochemical  potentials  of  electrons  at 
the  two  electrodes  in  terms  of  the  difference  in  chemical  po¬ 
tentials  of  oxygen  at  the  two  electrodes,  the  ionic  current  den¬ 
sity,  and  the  net  ionic  resistance,  is  applicable.  Note  that  there 
are  no  terms  containing  Ie  in  Eq.  (37).  However,  it  is  impor¬ 
tant  to  note  that  Ie  is  not  zero,  although  usually  \Ie  \  <<C  |/;|  in 
the  present  case  of  fuel  cells.  Thus,  the  net  measured  current 
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corresponding  to  f  +  Ie 

Fig.  6.  A  schematic  illustration  of  a  fuel  cell  showing  the  directions  of  oxy¬ 
gen  ion  and  electron  transport  through  the  membrane.  Electronic  and  ionic 
currents  are  in  the  opposite  directions. 


is  essentially  the  ionic  current.  Also,  note  that  f  <  0.  Thus, 
writing  I  as  the  magnitude  of  the  measured  current,  note  that 
Ii  ~  —I.  Then  Eq.  (37)  becomes 


II 


I{ipi  +  r'i  +  r"} 


(51) 


When  the  fuel  cell  is  at  open  circuit,  the  net  current  flow¬ 
ing  through  the  external  circuit,  7=0,  which  corresponds  to 
the  Nernst  voltage,  Eq,  given  by  Eq.  (39).  This  is  equiva¬ 
lent  to  the  establishment  of  (near)  global  equilibrium  for 
the  electrochemical  potential  of  oxygen  ions,  jl02- ,  that  is 
V£0 2-  0.  Note  of  course  that  V/z0 2-  is  not  identically 
zero,  and  therefore,  the  ionic  current  is  not  identically  zero 
either.  This  situation  is  now  equivalent  to  an  MIEC,  albeit 
with  a  very  low  electronic  transference  number,  and  thus 
Ie  =  — Ii ,  but  with  a  very  small  magnitude.  Now,  consider 
the  following  seven  cases  corresponding  to  the  open  cir¬ 


cuit  condition:  Case  E  re  r',  r"  Case  IE  r’  re,  r"  Case 


J! 
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e  ~  re.  The  vari- 


Case  VI:  re,  re  r",  and  Case  VIE  re  ~ 
ations  of  jl02- ,  cp  and  /xo2  for  the  seven  cases  are  shown  in 
Fig.  7.  All  seven  cases  correspond  to  the  open  circuit  con¬ 
dition,  with  [lpe  +  r'e  +  r"e]  ^  [Ipi  +  r\  +  r"],  and  therefore 
V£0 2-  %  0.  The  corresponding  difference  between  cpl  and 
<p11  for  all  seven  cases  is  nearly  equal  to  the  Nernst  voltage, 
given  by  Eq.  (39).  Note,  however,  that  the  details  of  the  vari¬ 
ations  of  pio2  and  <p  through  the  interfaces  and  membrane  de¬ 
pend  upon  the  relative  values  of  electronic  resistances  of  the 
interfaces  and  the  membrane.  From  the  standpoint  of  mem¬ 
brane  stability  in  a  low  po2  atmosphere  (stability  against 
reduction),  the  most  preferred  case  is  that  depicted  in  case 
III,  which  shows  that  if  r"  re,  r'e,  the  /jl o2  just  inside  the 
anode/electrolyte  interface  (that  is  /Xq2)  will  be  close  to  the 

cathode  side  chemical  potential  of  oxygen,  namely  /Xqo  .  This 
means,  under  such  conditions,  in  principle,  it  may  be  possible 


to  use  an  otherwise  unstable  electrolyte  in  fuel  atmosphere, 
such  as  rare  earth-oxide  doped  bismuth  oxide  for  example, 
for  a  fuel  cell  application. 1 1  All  cases  considered  correspond 
to  the  open  circuit  condition,  such  that  7/  +  Ie  =  0.  Thus,  this  is 
exactly  the  same  as  an  MIEC  membrane,  and  thus  po2  varies 
monotonically  from  the  value  at  one  electrode  to  the  other. 
When  the  fuel  cell  is  shorted: 

<pl  —  (pu  =  0  or  <pl  =  cp11 


Thus,  in  this  case,  the  electronic  current  density  is  zero,  i.e., 
Ie  =  0.  The  corresponding  short  circuit  current  density,  which 
is  exclusively  due  to  ionic  transport  through  the  membrane, 
is  given  by 


Po2  -  A  Pf  -  A 


4  e{L  +  !ff\  4e{£pi+r'i+r "} 

l  Sig"  / 


(52) 


Since  7/  <  0,  it  is  clear  from  Eqs.  (27)  and  (28)  that  in  this 
case  also  the  /xo2  varies  monotonically.  The  corresponding 
variations  of  jl02- ,  cp  and  \±o2  for  the  case  corresponding  to 
r’e  ~  re  ~  re  are  shown  in  Fig.  8.  Various  cases  can  be  readily 
examined  for  different  relative  values  of  electronic  specific 
resistances,  similar  to  Fig.  7.  It  can  be  readily  shown  that  for 
a  given  overall  performance  (voltage  versus  current  density), 
the  spatial  variations  of  po2  and  (p  can  vary  over  a  wide  range, 
depending  upon  the  relative  values  of  electronic  resistances. 

Now  let  us  examine  the  spatial  variations  of  /xo2  and  <p 
across  a  fuel  cell  membrane  when  the  fuel  cell  is  under  load. 
The  chemical  potentials  of  oxygen  just  inside  the  interfaces 
are  given  earlier.  Note  that  7/  <  0  but  Ie  >  0.  Thus,  the  term 
(r'jli  —  r'e I e)  can  be  replaced  by  —  (r- 1 7/ 1  +  r'e  \  Ie  |).  An  impor¬ 
tant  point  is  that  even  when  \Ie  \  |7/| ,  it  is  possible  for  r  ■  1 7/ 1 

and  r'e\Ie  \  to  be  of  comparable  magnitudes,  provided  <  »  5 
That  is,  the  contribution  of  Ie  for  the  estimation  of  /xo2  cannot 
be  a  priori  ignored.  Thus,  the  chemical  potentials  of  oxygen 
just  under  the  interfaces  are  given  by 


mo2  =  A  - 4e 


Se  I  h  I  +  g'i\Ii 


e 


SeS'i 


=  P02  ~  +  re\!e  I)  <  Po2 


(53) 


and 


(54) 


and  that  the  contribution  of  Ie  to  the  chemical  potential  of 
oxygen  cannot  be  a  priori  ignored.  Thus,  it  is  clear  that 
pil0?  >  /Xq9  >  /Xq9  >  /Xqo  .  That  is,  for  all  cases  (open  cir¬ 
cuit,  short  circuit,  and  under  load)  of  a  fuel  cell,  the  /xo2  varies 
monotonically  between  the  two  end  values.  Also,  therefore, 
as  long  as  /Xq2  (or  more  accurately  Mq2)  is  greater  than  the 
decomposition  chemical  potential  of  the  solid  electrolyte, 


1 1  There  are  several  difficulties,  however,  which  may  make  this  impractical. 
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Fig.  7.  Schematic  variations  of  jl02- ,  (p  and  pto 2  for  different  values  of  electronic  specific  resistances  at  open  circuit  through  a  fuel  cell.  Case  I:  re  »  r'e,  r'f 
Case  II:  r'eJp>  re,  r",  Case  III:  r"e  re,r'e,  Case  IV:  r'e ,  rne  ^  re,  Case  V:  re,  r"  r'e,  Case  VI:  re,r'e  r",  and  Case  VII:  r'  ~  re  ~  r".  For  all  seven  cases 
shown  above,  Re  >  >  /?;,  and  A/xo,  =  ^q2  —  ^o2  =  hxecf  which  leads  to  A (p  =  cpl  —  cpu  ~  hxed,  and  V/x02-  ~  0  (near  global  equilibrium  for  O2-). 
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Fig.  8.  Schematic  variations  of  jl02- ,  ip  and  ^cb  for  a  shorted  cell.  Note  that 
for  a  shorted  cell,  A (p  =  (p1  —  <p11  =  0,  and  the  electronic  current  through  the 
cell  is  zero. 

/4e2C°mp,  decomposition  of  the  electrolyte  cannot  occur.  It 

is  also  clear  that  cpl  >  cp' >  cp"  >  cp11 .  Schematic  variations  of 
the  chemical  potential  of  oxygen,  the  reduced  electrochem¬ 
ical  potential  of  electrons,  and  the  electrochemical  potential 
of  oxygen  ions  for  a  fuel  cell  under  load  are  shown  in  Fig.  9 
for  the  case  corresponding  to  r'e  ~  r”e  ~  re.  Note  that  this 
is  qualitatively  similar  to  that  for  the  case  of  MIEC  oxygen 
separation  membrane  shown  in  Fig.  3.  Similar  cases  can  be 
considered  for  other  values  of  electronic  area  specific  resis¬ 
tances. 

Similar  to  the  discussion  on  MIEC  oxygen  separation 
membranes,  the  corresponding  equivalent  circuit  for  a  fuel 
cell  is  given  in  Fig.  10,  with  once  again  the  provision  that  in 
steady  state  one  must  have  7/  =  /■  =  7"  and  Ie  =  /'  =  7". 
The  equivalent  circuit  is  similar  to  that  for  the  MIEC  mem¬ 
brane,  with  the  exception  that  there  is  now  an  external  load, 
whose  resistance  is  given  by  Rl  (in  ^  cm2).  The  internal 
EMFs,  and  their  relation  to  the  net  oxygen  chemical  poten¬ 
tial  difference  are  the  same  as  for  the  MIEC  membrane.  How¬ 
ever,  the  currents  are  different,  by  virtue  of  the  presence  of 
an  external  load.  It  is  easy  to  show  that 

E(Re  +  Rl) 

Ii  = - v  ; -  (55) 

RLRe  +  Ri(Re  +  RL) 


and 


ERe 

Rl  Re  +  Ri(Re  +  Rl) 


(57) 


In  the  above,  the  units  of  //,  Ie  and  II  are  in  A  cm-2 .  Note  that 
as  Rl  ->  oo,  7Z  and  Ie  are  identical  in  magnitude  (but  opposite 
in  sign,  of  course),  and  given  by  Eq.  (47),  which  is  the  same 
as  for  an  MIEC  membrane.  The  load  current,  II ,  is  given  by 


I L  =  h  +  Ie 


(58) 


The  overpotentials.  For  a  fuel  cell  operating  with  a  finite, 
nonzero  current  flowing  through  the  external  circuit,  there 
are  losses  associated  with  electrodes  as  well  as  the  ohmic 
loss  through  the  cell.  When  no  load  is  connected,  some  loss 
continues  to  occur  due  to  internal  electronic  leakage.12  The 
rate  of  external  work  done  or  power  is  given  by  l\  Rl  =  VI l 
where  V  is  the  voltage  across  the  load.  The  rates  of  losses 
through  the  electrodes  and  through  the  electrolyte  can  be 
similarly  calculated. 

For  example,  for  electrode  I  (cathode),  the  rate  of  loss 
(degradation  of  potential  work  into  heat — an  irreversible  pro¬ 
cess)  is  given  by 

fA  +  if’,  (59) 


If  r  •  and  r'e  are  not  constant,  but  depend  upon  the  local  chem¬ 
ical  potential  of  oxygen,  then  this  aspect  will  have  to  be 
included  in  the  analysis.  For  the  purposes  of  the  present  dis¬ 
cussion,  r\  and  r’e  are  assumed  to  be  constant.  The  usual  ap¬ 
proach  is  to  define  overpotential  at  an  electrode  as  a  measure 
of  the  loss  of  useful  voltage  into  irreversible  processes.  The 
implication  is  clearly  that  the  product  of  the  measured  overpo¬ 
tential  and  the  measured  current  (which  naturally  is  through 
the  external  circuit)  is  the  rate  of  loss  of  potential  work  into 
an  irreversible  process  (heat  generation).  Since  the  measure¬ 
ment  of  current,  namely  7^,  is  a  straight  forward  matter,  we 
will  define  overpotential  at  electrode  I  by  rjf  such  that 


rflL  =  Pr[  +  l}r[ 


e  e 


(60) 


ERl 

Rl  Re  +  Ri(Re  T-  Rl) 


(56) 


Fig.  9.  Schematic  variations  of  jl0 2- ,  ip  and  across  a  fuel  cell  under 
load. 


Thus: 


/V  1  /V 

^  _  1l  1  '  e '  e 


II 


fir(  +  7V 

11  1  e  e 


ii  +  /, 


e 


(61) 


In  the  above,  the  overpotential  at  electrode  I,  rjf ,  is  defined 
(not  necessarily  measurable)  in  terms  of  the  rate  of  loss  of 
potential  useful  work  and  the  measured  current  in  the  external 
circuit,  7^.  Since  the  rate  of  loss  of  potential  work  must  be 
equal  to  or  greater  than  zero  (a  measure  of  the  degree  of 
irreversibility),  the  product  ij'Il  must  be  positive.  In  terms 
of  the  measurable  parameters,  the  ionic  and  the  electronic 


12  In  a  typical  fuel  cell  with  negligible  electronic  leakage,  this  loss  is  neg¬ 
ligible,  yet  cannot  be  identically  zero.  However,  if  the  electrolyte  is  based 
on  ceria,  internal  electronic  leakage  can  be  significant. 
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!» 


Fig.  10.  An  equivalent  circuit  for  a  fuel  cell  under  an  external  load.  The  circuit  elements  between  the  two  filled  elliptical  symbols  are  not  physically  separable; 
they  collectively  represent  the  equivalent  circuit.  All  E’s  (Er ,  £;,  E")  are  of  the  same  sign. 


currents  are  given,  respectively,  by 


E  ILRe 
Ri  +  Re  Ri  +  Re 


and 


E 

Ri  +  Re 


+ 


hRi 

Ri  +  Re 


(62) 


(63) 


Note  that  as  Ie  — >  0  (open  circuit),  Ie  =  —f.  Substitution  for 
Ii  and  Ie,  respectively,  from  (62)  and  (63)  into  (61)  and  rear¬ 
rangement  gives 


T) 


! 


+  mRjf  ~  Ref)  Eff  +  Q 

(Ri  +  Re)2  (Ri  +  Re)2  IL(Ri  +  Re)2 

(64) 


As-defined,  the  parameters,  which  can  be  measured,  at  least 
in  principle,  are  the  ionic  and  the  electronic  area  specific 
resistances  of  the  electrolyte  and  interfacial  regions,  and  the 
current  flowing  through  the  external  circuit,  /l,  while  rj'  is  cal¬ 
culated  using  Eq.  (64).  Interesting  conclusions  emerge  from 
Eq.  (64).  Consider  first  the  short  circuit  limit.  For  this  case, 
Ie  =  013  and  Ii  =  —  Js  Substitution  into  Eq.  (60)  or  (64)  and 

simplification  gives  ^(short  circuit)  =  —  This  is  simply 
the  net  current  flowing  through  the  external  circuit  multiplied 
by  the  ionic  charge  transfer  resistance  of  electrode  I. 

Now  consider  the  open  circuit  case.  Under  open  circuit 
condition,  no  net  current  flows  through  the  external  circuit, 
that  is  Il  =  0.  As  no  current  flows  through  the  external  cir¬ 
cuit,  the  obvious  expectation  is  that  the  overpotential  is  zero. 


However,  as  discussed  in  what  follows,  this  conclusion  is 
erroneous.  This  is  because,  as  long  as  the  electronic  conduc¬ 
tivity  is  not  identically  zero,  it  means,  Ie  =  —Ii  ^  0.  Thus,  at 
open  circuit  condition,  according  to  the  definition  of  rj'  given 
inEqs.  (60),  (61)  and  (64),  rj'  — >  —  oo  or  \  rj'\  — >  oo.  This  looks 
like  a  strange  conclusion.  An  important  point  to  note,  how¬ 
ever,  is  that  the  relevant  physical  quantity  is  not  the  overpoten¬ 
tial,  rj' •  The  relevant  physical  quantity  is  v[Il\  which  is  the  rate 
of  loss  of  potential  useful  work  as  heat  (an  irreversible  pro¬ 
cess)  through  electrode  I  (cathode).  Note  that  while  rj'  ->  —  oo 

or  |  rf  I  — >  oo  as  Ie  — >  0,  the  product  rj'  I l  — >  —  g2 — is  a  fi- 

\Re5~Ri) 


nite  positive  limit — as  >  0.  But  this  is  simply  the  rate  of 

loss  of  potential  work  as  heat  through  electrode  I  due  to  inter¬ 
nal  leakage.  This  is  the  significance  of  overpotential  at  open 
circuit. 

An  examination  of  Eq.  (64)  shows  that  the  dependence  of 
rj'  on  I i  is  nonlinear.  This  is  despite  the  fact  that  all  processes 
are  assumed  to  be  describable  by  various  resistances.14  The 
limit  of  Eq.  (64)  as  Re  ->  oo  gives15 


(65) 


which  is  the  usual  result.  But  this  is  valid  only  if  the  elec¬ 
tronic  resistance  is  infinite  (very  large  compared  to  the  ionic 
resistance).  Also,  only  in  this  case  the  rj'  is  the  same  as  the 
experimentally  measured  overpotential. 

The  preceding  suggests  that  a  plot  of  Ie  versus  rj'16  should 
strictly  not  go  through  (0,  0)  coordinates.  Eq.  (64)  shows  that 
there  should  be  a  minimum  in  \rj'\.  This  minimum  occurs  at 


13  That  is,  the  electronic  current  is  identically  equal  zero.  This,  however, 
does  not  imply  that  local  equilibrium  is  violated.  The  reason  the  electronic 
current  is  zero  in  this  case  is  because  the  electrodes  are  shorted.  In  such  a 
case,  even  though  Ie  through  the  membrane  is  zero,  that  through  the  combined 
system  comprising  the  membrane  and  the  external  short  is  not  zero.  Also, 
the  preceding  is  only  for  steady  state.  If  the  steady  state  has  not  been  estab¬ 
lished,  the  position  and  time-dependent  electronic  current,  Ie(x,  t ),  through 
the  membrane  need  not  be  zero,  even  when  externally  shorted. 


14  Note,  however,  that  the  relationship  between  voltage  across  the  cell 

F  R 

and  the  measured  current  is  still  linear  and  is  given  by  V  =  {R+R  )  — 

(D  .  D  \ 

R  l+R  \  \Ii\.  \fRe  oo,  this  reduces  to  the  usual  equation  V=E  —  Ri\Ii\. 
15  Once  again,  to  emphasize,  Re  strictly  cannot  go  to  infinity  since  this 
violates  the  assumption  of  local  equilibrium.  The  above  simply  implies  that 
Re  is  very  large. 

16  To  reemphasize,  the  as-defined  if  is  calculated,  not  measured. 
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a  current  density  given  by 


| /l  (min)  |  =  E 


(66) 


Note  that  when  Re  — >►  oo,  |/zXmin)|  — >►  0.  Note  also  that  al¬ 
though  |  rj'\  exhibits  a  minimum,  the  product  I  if  monoton- 
ically  increases  with  increasing  \Il\,  from  its  lowest  value 
at  OCV  to  its  highest  value  at  short  circuit.  This  means,  the 
rate  of  loss  of  potential  useful  work  due  to  irreversibility  of 
the  process  associated  with  overpotential,  monotonically  in¬ 
creases  with  load  current. 

For  electrode  II  (anode)  and  the  bulk  membrane,  similar 
equations  can  be  written.  They  are  as  follows: 


For  electrode  II  (anode): 


„  =  (J&l  +  R}QIL  2 E(Rjrf:  -  Rer ") 

(Ri  +  Re)1  ( Ri  +  Re)2 

E2(r"  +  r") 

+  IL(Ri  +  Re)2 


In  the  short  circuit  limit 


0  (short  circuit)  = 


and  under  open  circuit  conditions 


E2(r"  +  r"e) 

(Re  +  Rif 


as 


0 


The  limit  of  Eq.  (67)  as  Re  — >  oo  gives  r f  =  /^r". 
For  the  bulk  membrane: 


(Rfn  +  Rjre)iL  2E(Rire  -  Ren) 
(Ri  +  Re)2  (Ri  +  Re)2 

lL(Ri  +  Re)2 


(67) 


(68) 


(69) 


Note  that  when  Re  — >►  oo,  rj  — >►  IlRu  From  Eq.  (72),  the  rate 
of  loss  of  potential  work  into  heat  (irreversible  process)  in  all 
parts  of  the  cell  is  given  by 


ReRi  2 
— — /?  + 

Ri  +  Re 


E 2 

Ri  +  Re 


(73) 


which  for  Re  ->  oo  reduces  to  l\  Ri . 

Overpotentials  in  terms  of  thermodynamic  potentials.  The 
overpotential  can  be  expressed  in  terms  of  chemical  potential 
of  O2,  /xo2,  and  reduced  (negative)  electrochemical  potential 
of  electrons,  ip,  by  substituting  for  f  and  Ie  into  Eq.  (60),  and 
using  Eqs.  (29)  and  (30).  In  general,  rather  complicated  equa¬ 
tions  result.  In  what  follows,  only  the  limiting  case  wherein 
the  electronic  resistance  is  much  greater  than  the  ionic  re¬ 
sistance  will  be  examined,  that  is,  the  case  corresponding  to 
Re  — >  00.  For  this  case,  it  can  be  readily  shown  that  overpo¬ 
tential  at  electrode  I  is  given  by 


0  = 


r,H  _  r,  I 

Mq2-  d02- 

2e 


r- 

1 

R[ 


(74) 


given  in  terms  of  the  net  drop  in  electrochemical  potential 
of  oxygen  ions  across  the  cell,  the  ionic  resistance  of  the 
interface  at  I,  and  the  net  ionic  resistance,  or 


-  Aq 2- 
2e 


(75) 


given  in  terms  of  the  drop  in  electrochemical  potential  of 
oxygen  across  electrode  I.  As  defined,  the  f  is  negative,  and 
the  net  current,  /l,  is  also  negative.  Similar  equations  can  be 
written  for  f  and  rji  in  terms  of  electrochemical  potentials  of 
oxygen  ions.  Note,  however,  that  Eqs.  (74)  and  (75)  are  valid 
only  if  the  net  electronic  resistance  is  much  larger  than  the 
net  ionic  resistance.  Also,  only  in  such  a  case  the  measured 
overpotential  is  the  same  as  v[  defined  here. 


2.4.  Oxygen  separation  under  the  application  of  an 
external  voltage 


In  the  short  circuit  limit 


ip  (short  circuit)  = - r; 

Ri 


and  under  open  circuit  conditions 


mh 


Epin  +  re) 

(Re  +  Ri)2 


as 


(70) 


The  limit  of  Eq.  (69)  as  Re  — >  00  gives  r)  =  ////. 

The  total  polarization  loss  is  the  sum  of  the  three,  namely 


rj  =  v[  +  rji  +  0 

which  is  also  given  by 

Re  R,  E2 

n  =  — - — - — 1L  4 - 

Ri  +  Re  E(R,  +  Re) 


(71) 


(72) 


Such  would  be  the  case  with  yttria- stabilized  zirconia 
(YSZ)  or  rare  earth  oxide-doped  ceria  as  the  electrolyte,  with 
suitable  electrodes  applied  on  the  two  surfaces.  In  the  first  two 
cases  considered,  namely  an  MIEC  membrane  and  a  fuel  cell, 
the  A/xo2  and  are  related  to  each  other.  Specifically,  upon 
the  application  of  A/xo2  across  the  membrane,  a  correspond¬ 
ing  A cp  appears  across  the  membrane.  In  the  present  case,  an 
external  voltage  is  applied  across  the  two  electrodes,  which 
effectively  fixes  the  A  ip  across  the  membrane.  This  provides 
an  additional  experimental  degree  of  freedom.  Fig.  1 1  shows 
a  schematic  illustration. 

In  the  present  situation,  note  that  oxygen  is  being  pumped 
through  the  membrane  from  left  to  right.  Thus,  f  <  0.  Simi¬ 
larly,  electrons  are  also  being  pumped  through  the  membrane 
from  left  to  right.  Thus,  Ie  <  0.  That  means,  both  f  and  Ie  are 
of  the  same  sign,  unlike  the  previous  two  cases.  In  this  case, 
both  A/xq2  and  A  ip  are  experimental  parameters,  which  can 
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corresponding  to  /,•  +  Ie 


Fig.  1 1.  A  schematic  illustration  of  an  oxygen  separation  membrane  under 
an  externally  applied  voltage.  Electronic  and  ionic  currents  are  in  the  same 
direction. 


be  arbitrarily  selected.  The  general  transport  equations  are 
the  same  as  before,  and  the  net  electronic  and  ionic  current 
densities  are  given  by  Eqs.  (29)  and  (30),  respectively. 

The  oxygen  ion  current,  as  before,  is  given  by  Eq.  (23), 
and  the  chemical  potentials  of  oxygen  just  under  the  two 
interfaces  are  given,  as  before,  by  Eqs.  (27)  and  (28).  Since 
now  both  li  <  0  and  Ie  <  0,  Eqs.  (27)  and  (28)  can  be  written 
as  follows: 


(76) 


(77) 


Note  the  negative  sign  in  parentheses.  Thus,  if  r 

note  that  fil0i  <  /Xq2  ,  and  if  r •  |  /*■  |  <  r'e  \  Ie  | ,  note  that  /Xq2  > 

/Xqo  .  Similarly,  if  r"  |7)|  >  r"e  \  Ie  \ ,  note  that  /Xqo  >  /Xq2  ,  and 

if  r’l  |  f  |  <  r"e  1 7e  | ,  note  that  /1q2  <  p1^. 

In  the  following  discussion,  it  will  first  be  assumed  that 
r-\Ii\  >  r'e\Ie\  and  r-'|/;|  »  r"e\Ie\.  Thus,  note  that  /Xq2  < 

/Xqo  and  /Xqo  >  /Xq0.  Substitution  of  Eq.  (30)  into  Eq.  (76) 
gives 


Mo,  =  Mo, 


(Mq2  ~  PfVj  Ae{<pl  -  <?n}r- 
ipi  +  r  •  +  r"  tpi  +  r(-  +  r" 


(78) 


Suppose  the  chemical  potential  of  oxygen  is  the  same  at  both 
electrodes,  that  is  /Xq2  =  /Xq2  .  Then: 


4  e{(pl  —  <pll}r[ 
Ipi  +  r[  +  r’! 


4  e{cp11  —  (pl}r't 
ipi  +  r\  +  r’! 


(79) 


Since  ^n>^x,  note  that  p!0n  <  /Xq2.  The  absolute  value  of 
/Xqo  in  relation  to  /Zqo  depends  upon  the  various  transport  pa¬ 


rameters,  especially  r  •,  and  the  net  voltage  difference,  namely, 
A  cp  =  (p11  —  (pl.  The  higher  are  the  r-  and  Acp  =  cp11  —  <pl,  the 
lower  is  .  If  /Xq?  is  lower  than  the  decomposition  poten¬ 
tial  of  the  electrolyte,  decomposition  can  occur,  effectively 
increasing  the  magnitude  of  the  ionic  current  density,  beyond 
that  given  by  Eq.  (23).  In  the  case  of  zirconia,  for  example, 
the  decomposition  reaction  may  involve  the  following  steps. 

When  the  p!0i  <  /XQ2comp,  reaction: 

ZrC>2  — >  Zr4+  +  202- 

occurs  just  under  electrode  I.  The  oxygen  ions  released  trans¬ 
port  through  the  membrane  to  electrode  II  where  the  follow¬ 
ing  reaction  occurs: 

202“  02(g)  +  4e' 

The  electrons  released  transport  in  the  external  circuit  to  elec¬ 
trode  I  where  the  following  reaction  occurs: 

Zr4+  +  4e  — >  Zr(metal) 

The  net  reaction  is 


ZrC>2  — >  Zr(metal,  formed  at  negative  electrode) 
+  02(g,  released  at  positive  electrode) 


The  corresponding  current  appears  in  the  external  circuit,  in 
addition  to  the  current  due  to  oxygen  pumping  from  left  to 
right.  An  important  point  is  that  the  corresponding  electron 
transport  occurs  in  the  external  circuit. 

Also  an  important  point  to  note  is  that  only  in  a  very  special 
case  would  one  have  /Zq2  =  Mq?.17  In  almost  all  practical 
situations,  such  as  the  one  assumed  here,  one  would  almost 
always  have  /i'0i  <  p}0  (for  the  assumed  case  of  r'i\U\  » 
r'e\Ie\). 

Similar  analysis  can  be  carried  out  for  the  interface  at  II. 
For  the  case  of  /Xq2  =  /Xq2  ,  it  is  seen  that 


4  e{cp11  —  (pl}r'j' 
ipi  +  r[  +  r" 


(80) 


for  the  assumed  case  of  r 1  f  \  r"e\Ie\.  That  is,  one  would 
have 


Mo2  >  mS2  =  Mo2 

Schematic  variations  of  /xo2  and  cp  for  this  special  case  are 
shown  in  Fig.  12.  Also  shown  in  the  figure  is  a  schematic 
variation  of  jl02- . 

If,  on  the  other  hand,  r  •  7/ 1  <  r'e  \  Ie  \  and  r'(  \  U  \  <  r"e  \  Ie  \ ,  we 
would  have  had  /Xqo  >  /Xqo  and  /Xqo  <  /Xqo  .  In  such  a  case, 
the  variation  of  fi o2  and  (p  would  be  as  shown  schematically  in 
Fig.  13.  In  such  a  case,  the  decomposition  of  the  electrolyte 
can  occur  at  the  positive  electrode  instead.  Once  again  the 


17  For  the  general  case  (Eq.  (76)),  for  /Xq2  to  be  equal  to  /j1^  ,  it  is  necessary 
that  P |  f  |  —r'e\Ie\  =0.  Clearly,  this  is  a  very  rare  case,  and  unlikely  to  be 
realized  in  practice. 
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Fig.  12.  Schematic  variations  of  /xo9 ,  <p,  and  jl0 2-  across  an  oxygen  separa¬ 
tion  membrane  under  an  externally  applied  voltage.  It  is  assumed  here  that 
r't  \f  \  >  r'e\Ie  \  and  r" \f  \  >  r" \Ie\.  Note  that  the  /x02-  decreases  from  left  to 
right,  consistent  with  net  oxygen  transport  from  left  to  right. 

reaction  steps  are  as  follows.  When  the  /Xqo  <  /XQ2comp,  re¬ 
action 

ZrC>2  — >  Zr4+  +  202- 


Fig.  14.  Schematic  variations  of  pto2,  cp,  and  jl0 2-  across  an  oxygen  sep¬ 
aration  membrane  under  an  externally  applied  voltage.  It  is  assumed  here 
that  m  >  r'e\Ie  \  and  r"\li\  <  <  141  .  In  this  case,  the  \iq2  in  the  membrane 
everywhere  is  lower  than  that  in  the  gas  phase.  Whether  the  /io2  increases 
or  decreases  from  left  to  right  depends  upon  relative  magnitudes  of  transport 
parameters.  Note  that  the  jl02-  decreases  from  left  to  right,  consistent  with 
net  oxygen  transport  from  left  to  right. 


occurs  just  under  electrode  II.  Also  at  electrode  II 
202-  -*  02(g)  +  4e' 

Electrons  released  transport  through  the  external  circuit  to 
electrode  I,  and  then  through  the  membrane  towards  electrode 
II,  where  the  following  reaction  occurs: 

Zr4+  +  4er  — >  Zr(metal) 

The  net  reaction  is 

ZrC>2  — >  Zr(metal,  formed  at  positive  electrode) 

+  C>2(g,  released  at  positive  electrode) 

In  this  case,  however,  it  is  necessary  for  the  corresponding 
electron  transport  to  occur  through  the  membrane.  Since  the 
membrane  is  predominantly  an  ionic  conductor,  electronic 
conduction  is  expected  to  be  negligible.  The  implication  is 
that  although  the  analysis  shows  that  decomposition  can  oc¬ 
cur  at  the  positive  electrode,  the  associated  transport  kinetics 


Fig.  13.  Schematic  variations  of  /xo9 ,  <p,  and  jl0 2-  across  an  oxygen  separa¬ 
tion  membrane  under  an  externally  applied  voltage.  It  is  assumed  here  that 
P \Ii\  <  r'e\Ie  \  and  r'! \f  \  <  r"e  \ Ie | .  Note  that  the  /xQ 2-  decreases  from  left  to 
right,  consistent  with  net  oxygen  transport  from  left  to  right. 


may  be  very  sluggish.  Thus,  in  practice,  this  case  may  not  be 
observed. 

Additional  cases  can  be  similarly  considered.  They  are  as 
follows: 

Suppose  r[\Ii\  <  rfe\Ie\,  so  that  \i'0i  >  /Xq2,  and  r"|//l  < 

r "  1 7e  | ,  so  that  /Xqo  <  /Xqo  .  Fig.  14  shows  the  expected  vari¬ 
ation  of  the  chemical  potential  of  oxygen  through  the  mem¬ 
brane  for  /Xq2  =  .  Note  that  in  this  case,  the  fio2  in  the 

membrane  is  lower  than  that  in  the  gas  phase. 

Suppose  now  r •  |  f  \  >  r'e\ I e\,  so  that  /Hq2  <  /Xq2,  and 

rf{ |  f  |  >  r”\Ie\,  so  that  /Xq2  >  /Xq9.  Fig.  15  shows  the  ex¬ 
pected  variation  of  the  chemical  potential  of  oxygen  through 
the  membrane.  In  this  case,  the  /xo2  in  the  membrane  is 
higher  than  that  in  the  gas  phase.  The  preceding  shows  the 
importance  of  incorporating  the  electronic  current  through 


Fig.  15.  Schematic  variations  of  /xcb,  (p ,  and  jl0 2-  across  an  oxygen  sep¬ 
aration  membrane  under  an  externally  applied  voltage.  It  is  assumed  here 
that  m<r'e\  Ie  \  and  >  re  I  h I  .  In  this  case,  the  pto2  in  the  membrane 
everywhere  is  higher  than  that  in  the  gas  phase.  Whether  the  /xo0  increases 
or  decreases  from  left  to  right  depends  upon  relative  magnitudes  of  transport 
parameters.  Note  that  the  jl02-  decreases  from  left  to  right,  consistent  with 
net  oxygen  transport  from  left  to  right. 
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Fig.  16.  An  equivalent  circuit  for  an  oxygen  separation  membrane  under  an  applied  DC  voltage.  The  circuit  elements  between  the  two  filled  elliptical  symbols 
are  not  physically  separable;  they  collectively  represent  the  equivalent  circuit.  The  E’s  ( E ' ,  Ei ,  E")  can  be  of  different  signs,  and  depend  upon  the  applied 
voltage,  Ea,  and  the  logarithm  of  the  ratio  of  oxygen  partial  pressures  at  the  two  electrodes.  Their  sum,  however,  is  independent  E\,  but  depends  only  on  the 
logarithm  of  the  ratio  of  oxygen  pressures  at  the  two  electrodes. 


bulk  and  also  across  interfaces,  however  small,  into  the 
transport  equations,  even  when  dealing  with  membranes 
made  of  predominantly  oxygen  ion  conductors. 


The  equivalent  circuit  for  an  oxygen  separation  membrane 
under  an  externally  applied  voltage  is  given  in  Fig.  1 6  with  the 
provision  that  in  steady  state  7Z  =  /•  =  7"  and  Ie  =  Ee  =  7". 
The  same  equivalent  circuit  is  applicable  for  the  transient 
case,  in  which  the  ionic  and  electronic  currents  are  functions 
of  position  and  time,  but  Kirchoff’s  equations  continue  to 
hold.18  In  what  follows,  only  the  steady  state  is  considered. 
In  this  case,  as  before: 


E  =  E’  +  Ei  +  E” 


RT 

H - In 

4  F 


+ 


RT 
4 ~F 


In 


RT .  (po2\  kBr  (p^A 

^ln  U  j = ^ln  UJ 

(81) 


internal  EMFs  at  steady  state  are  given  by 


and 


(82) 


Note  that  the  terms  in  parentheses  determine  the  signs  of  the 
individual  EMFs,  yet  the  sum  of  the  three  EMFs  is  zero. 

If  /Xq2  ^  /Xq2,  then  E^  0.  Then,  it  can  be  readily  shown 
that  the  various  EMFs  are  given  by 


Ei  =  Ea 
E”  =  Ea 


+  E 


Ri 


+  E 


Ri 


+  E 


r± 

Ri 


and 


(83) 


In  the  previous  two  examples,  that  is  an  MIEC  membrane 
and  a  fuel  cell,  all  three  internal  EMFs  were  of  the  same 
sign — consistent  with  a  monotonic  variation  of  /jLo2  from  one 
side  to  the  other.  In  the  present  example  of  voltage-driven 
oxygen  separation,  however,  this  is  not  the  case.  Thus,  the 
signs  of  the  three  EMFs  in  general  are  different  as  discussed 
earlier.  An  important  point  to  note  is  that  the  net  E  in  the 
above  equation,  however,  is  still  independent  of  the  exter¬ 
nally  applied  voltage,  and  is  only  a  function  of  A/xo2  existing 
across  the  membrane.  However,  the  individual  internal  EMFs, 
namely,  E\  Ei ,  and  E"  are  functions  of  the  applied  voltage, 
Ea ,  and  thus  do  determine  the  stability  of  the  membrane. 
For  the  case  where  /Xq2  =  /Xq2  ,  the  E  =  0.  The  corresponding 


18  For  membranes  with  spatially  dependent  transport  properties  (but  with¬ 
out  the  incorporation  of  interfacial  effects),  the  transient  case  was  previously 
examined  in  [15,16]. 


and  the  sum  of  the  three  EMFs  equals  that  given  by  Eq.  (81). 

The  overpotentials.  In  the  case  of  oxygen  pumping  under 
an  applied  voltage,  an  external  agency  does  work  on  the  sys¬ 
tem.  Under  isothermal  conditions,  with  oxygen  gas  assumed 
to  behave  as  an  ideal  gas,  the  internal  energy  is  a  function 
of  temperature  only.  Thus,  no  change  in  the  internal  energy 
occurs  upon  isothermal  pumping.  We  will  first  consider  here 
the  case  wherein  the  oxygen  partial  pressures  at  the  two  elec¬ 
trodes  are  identical  so  that  /Xq2  =  /Xq2  and  the  net  E  given 
by  Eq.  (81)  is  zero.  This  means  all  work  done  on  the  system 
is  degraded  as  heat  (an  irreversible  process).  The  net  ionic 
and  electronic  current  densities  are  given  by  7Z  =  —  ^  and 

Ie  =  —  "I4.  Then,  the  rate  of  loss  of  input  work  into  heat  at 
electrode  I  is  given  by 


-1-  T^rr  —  F ^ 
E  *7  '  *e^e  ^ A 


(84) 
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Once  again,  the  overpotential  v[  is  defined  in  terms  of  the 
current  measured  in  the  external  circuit,  /,  by 


/?r-  +  /V 

i  i  1  e  e 


(85) 


where  I  =  —  Ea 


and  the  overpotential  is  given  by 


T) 


/ 


r  '■  R 2  +  r'  7?? 
i  e  1  e  i 

RiReiRi  +  Re) 


(86) 


In  this  case,  as  f  and  /g  are  of  the  same  sign,  the  rj '  is  the  same 
as  the  measured  overpotential.  If  the  electronic  resistance 
is  much  larger  than  the  ionic  resistance,  that  is  in  the  limit 
Re^  oo,  the  above  reduces  to 


ea  ,  , 

- ri  =  7r, 

R;  1  1 


(87) 


Working 

Electrode 


Reference 

Electrode 


Fig.  17.  A  schematic  illustration  of  the  three-electrode  system.  A  DC  voltage 
is  applied  across  the  working  electrode  and  the  counter  electrode  (on  the  op¬ 
posite  face).  Both  current  and  voltage  between  the  working  and  the  reference 
electrodes  are  measured  as  a  function  of  the  applied  voltage.  Usually,  the 
testing  is  conducted  in  air,  although  can  be  conducted  in  other  environments. 


In  terms  of  electrochemical  potential  of  oxygen  ions,  the  over- 
potential  can  be  given  as  before  by 


0 


! 


Aq2- 


2e 


(88) 


In  the  above,  as-defined,  both  I  and  rj'  are  negative.  Similar 
equations  can  be  given  for  77  and  77",  which  are,  respectively, 
overpotential  losses  across  the  bulk  membrane  and  electrode 
II.  It  is  readily  shown  that  the  net  rate  of  degradation  of  input 

work  into  heat  is  given  by  (rf  +  r 7  +  rj")I  =  E\  , 

e 2 

which  for  Re—*o o  reduces  to  . 

Suppose  we  now  consider  the  case  wherein  /Xq  7^  /r  J*,2 , 

such  that  E  =  — Ee — 1  7^  0.  In  such  a  case,  it  can  readily 
shown  that  the  ionic  and  electronic  currents  are  given  by 


Ea  -  E 

Ri 


and  /(. 


Ea 

Re 


(89) 


When  I  mol  of  oxygen  gas  is  pumped  from  pl0  to  plX  ,  the 
amount  of  potential  work  degraded  as  heat  is 


Q  =  4  F(Ea 


E)  + 


4 FEE A  R, 
Ea  -  E  Re 


(90) 


and  the  net  work  of  compressing  1  mol  of  gas  from  pf 

(Po2)  to  (Pq2)  is  §iven  by  w  =  4 FE  =  RT  In  • 

Note  that  if  Re  Ri  (or  Re  -a  oo),  the  above  reduces  to 
Q  ~  4F(Ea  —  E).  If,  however,  Re^o o,  but  if  E\  =  E ,  then 
the  Q  — >►  00.  Note  that  in  this  case,  =  0.  Thus,  it  will  take 
infinite  time  to  transport  1  mol  of  oxygen.  However,  Ie  ^  0. 
As  a  result,  infinite  electronic  charge  will  have  to  be  trans¬ 
ported  in  the  time  required  to  transport  1  mol  of  oxygen  gas, 
which  leads  to  Q  — >►  oo.  Finally,  the  overpotentials  can  be  eas¬ 
ily  determined  in  terms  of  Ie  and  the  various  area  specific 
resistances,  using  the  approach  described  earlier. 


2.5.  The  analysis  of  three -electrode  system 

Often,  the  placement  of  a  reference  electrode  in  devices 
such  as  solid  oxide  fuel  cells  (SOFC)  raises  questions  con¬ 
cerning  the  validity  of  the  measurement,  because  the  refer¬ 
ence  electrode  location  is  on  the  surface.  In  thin,  electrolyte 
film,  electrode- supported  SOFC,  this  leads  to  substantial  er¬ 
rors.  In  an  attempt  to  circumvent  these  uncertainties,  several 
investigators  have  used  the  three-electrode  system.  It  consists 
of  a  thick  electrolyte  disc,  upon  which  working  and  counter 
electrodes  are  placed  on  the  opposite  sides  of  the  disc  in  a 
symmetric  arrangement,  and  a  reference  electrode  is  placed  in 
a  symmetric  arrangement  with  respect  to  the  two  electrodes. 
Fig.  17  shows  a  typical  arrangement.  Then,  in  order  to  inves¬ 
tigate  electrode  kinetics  at  the  working  electrode,  which  may 
be  a  prospective  cathode  for  SOFC,  a  DC  voltage  is  applied 
across  the  working  and  the  counter  electrodes. 19  Both  current 
through  the  cell  and  voltage  across  the  working  and  reference 
electrodes  are  measured  as  functions  of  the  applied  voltage. 
Corrected  for  the  ohmic  loss,  the  measured  voltage  between 
the  working  and  the  reference  electrodes  gives  a  measure  of 
the  electrode  overpotential  at  the  working  electrode.  In  what 
follows,  implications  of  such  measurements  are  described 
using  the  analysis  presented  so  far. 

The  membrane  or  the  disc  is  predominantly  an  oxygen 
ion  conductor.  Thus,  note  that  \Ie  \  |C|.  In  the  following 

illustration,  it  is  assumed  that  r •  |  f  \  ^>>  r’e\Ie\  and  r'l\h I  » 


19  In  many  experimental  procedures,  the  instrument  is  operated  in  such  a 
mode  so  as  to  only  track  the  voltage  between  the  working  and  the  refer¬ 
ence  electrodes,  and  little  attention  is  paid  to  the  actual  voltage  between  the 
working  and  the  counter  electrodes.  It  is  to  be  emphasized  that  from  the 
standpoint  of  electrolyte  stability  and  the  present  discussion,  it  is  important 
to  know  the  actual  applied  voltage,  which  is  between  the  working  and  the 
counter  electrodes.  That  is,  the  measurement  of  voltage  between  the  working 
and  reference  electrodes  is  not  sufficient  to  describe  the  problem  completely. 

20  Since  \h  \  \Ie\,  this  also  implies  thatf^ 
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The  chemical  potential  of  oxygen  just  under  interface  I  is 
given  by  Eq.  (78).  The  testing  is  typically  done  in  air.  Thus, 
the  chemical  potential  of  oxygen  at  the  two  electrodes  (just 
outside  the  electrodes)  is  the  same.  If  concentration  polar¬ 
ization  is  negligible,  the  chemical  potentials  of  oxygen  at  the 
two  electrode/electrolyte  interfaces,  namely  /Xqo  and  /Xq0  are 

identical,  that  is  /Xq2  =  /Xqo  .  Let  us  consider  this  case  first. 
Eq.  (78)  then  becomes 


Mo2  =  Po2 


4  e{cp11  —  (pl}r'j 
Ipi  +  r  •  +  r" 


(91) 


If  /Xqo  decreases  below  the  decomposition  chemical  poten¬ 
tial,  /XQ2comp,  the  decomposition  can  occur  just  under  elec¬ 
trode  I,  which  manifests  as  excess  current.  It  is  also  possi¬ 
ble  that  a  nonstoichiometry  (oxygen  deficiency)  may  develop 
at  a  lower  applied  voltage,  which  can  also  lead  to  excess 
current  (beyond  that  due  to  the  cathodic  reaction,  namely 
^C>2  +  2e  — >*  O  ).  Since  the  experimenter  has  no  way  of 
determining  how  much  of  the  current  is  due  to  decomposition 
(or  nonstoichiometry  development)  and  how  much  is  actually 
attributable  to  the  cathodic  reaction  occurring  at  electrode  I, 
the  entire  current  is  inadvertently  assigned  to  the  cathodic 
reaction  leading  to  an  overestimation  of  cathodic  activity, 

and  thus  an  incorrect  result.  The  ^^comp  can  in  principle  be 
estimated  from  the  known  thermodynamic  data.  The  usual 
assumption  made  during  such  tests  is  that  all  that  needs  to  be 
done  is  to  correct  for  the  ohmic  drop,  in  order  to  ensure  that 
the  applied  voltage  is  less  than  some  critical  value  to  ensure 
that  the  chemical  potential  in  the  electrolyte  is  greater  than 

4C°mp  [21].  For  the  condition  /Xq2  ~  Eq.  (91)  can 

be  rearranged  as 


(0  —  w  ~ 


decomp 
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4e 


tPi  +  rt  +  r  ■ 
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(92) 


Eq.  (92)  shows  that  if  the  applied  voltage  exceeds  the  above 
value,  decomposition  of  the  electrolyte  will  occur  at  electrode 
I.  The  equation,  however,  shows  that  there  is  no  way  to  correct 
for  the  ohmic  loss  as  usually  both  r\  and  r"  are  unknowns. 
In  the  extreme  case  in  which  r\  r"  and  r\  Ipt,  Eq.  (92) 
reduces  to 


decomp 

P02 


4e 


(93) 


Eq.  (93)  gives  the  minimum  applied  voltage  necessary  to 
cause  decomposition  of  the  electrolyte  for  the  relative  values 
of  transport  parameters  assumed  here.  Beyond  this  voltage, 
the  measured  current  will  have  two  components;  current  due 
to  the  cathodic  reaction  and  current  due  to  decomposition. 
In  many  cases,  nonstoichiometry  may  appear  at  applied  volt¬ 
ages  lower  than  given  by  Eq.  (92)  or  (93),  and  this  also  will 
contribute  to  excess  current.  In  either  case,  the  measured  cur¬ 
rent  will  be  greater  than  that  corresponding  to  the  cathodic 


Fig.  18.  A  schematic  illustration  of  the  measured  current  (or  current  density) 
( — )  vs.  measured  overpotential  (corrected  for  the  ohmic  drop)  using  a  three- 
electrode  system  under  an  applied  DC  voltage.  It  contains  two  parts:  (a)  the 
cathodic  current  ( — ),  namely  due  to  the  reaction  ^  O2  +  2e'  — ►  02~ ,  and  (b) 
the  current  due  to  the  development  of  nonstoichiometry  or  decomposition 
( — ).  In  practice,  it  is  quite  possible  that  the  nonstoichiometry/decomposition 
current  may  be  much  greater  than  the  cathodic  current.  Under  such  condi¬ 
tions,  the  three-electrode  system  can  grossly  overestimate  the  cathodic  ac¬ 
tivity.  The  schematic  shows  both  current  and  overpotential  plotted  on  linear 
scales.  If  the  current  is  plotted  on  a  logarithmic  scale,  the  curves  will  be 
convex  up. 

reaction,  leading  to  an  overestimation  of  the  electrode  ac¬ 
tivity.  That  is,  one  may  erroneously  conclude  that  the  elec- 
trocatalytic  activity  of  the  prospective  cathode  is  far  better 
than  it  really  is.  Fig.  18  shows  a  schematic  illustration  of  the 
experimentally  measured  current  (or  current  density)  plot¬ 
ted  versus  experimentally  measured  overpotential  in  a  typi¬ 
cal  three-electrode  system  under  an  applied  DC  bias.  In  the 
schematic  illustration  shown  in  Fig.  18,  both  the  current  and 
overpotential  are  plotted  on  linear  scales.  Hence  the  curvature 
is  concave  up.  If  the  overpotential  is  plotted  on  a  linear  scale 
but  the  current  is  plotted  on  a  logarithmic  scale,  the  curvature 
will  be  convex  up. 

In  the  preceding  discussion,  it  was  assumed  that  concen¬ 
tration  polarization  is  negligible.  In  practice,  if  the  electrode 
thickness  is  large  and/or  electrode  porosity  is  low,  such  an 
assumption  may  not  be  justifiable.  A  review  of  literature 
shows  that  rarely  is  any  attention  paid  to  this  aspect  (thick¬ 
ness  and  porosity)  of  the  working  electrode  when  making 
such  measurements.  Suppose,  for  example,  that  concentra¬ 
tion  polarization  is  not  small  at  the  working  electrode,  cath¬ 
ode  (electrode  I).  In  such  a  case  under  an  applied  voltage, 
the  variation  of  oxygen  partial  pressure  through  electrode  I 
will  schematically  appear  as  shown  in  Fig.  19.21  In  such  a 
case,  /Xq?  <  /Xq2  where,  as  before,  /Xq2  is  the  chemical  po¬ 
tential  of  oxygen  just  outside  the  interface  between  electrode 
I  and  the  electrolyte.  In  the  extreme  case,  if  the  porosity  is  low 
and/or  the  electrode  thickness  is  large,  the  net  cathodic  current 
is  gas  phase  diffusion-limited  or  concentration  polarization- 
limited  (through  the  porous  electrode).  In  such  a  case,  the 
plQi  reaches  some  low  value,  such  that  beyond  some  applied 


21  Although  in  Fig.  19  the  variation  in  po2  is  shown  linear,  this  need  not 
be  the  case.  If  there  is  a  gradation  in  porosity,  the  variation  in  po0  will  be 
nonlinear. 
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Fig.  19.  A  schematic  variation  of  oxygen  chemical  potential,  /z<>,  (or  oxy¬ 
gen  partial  pressure,  po2 )  across  a  membrane  under  an  applied  DC  voltage. 
It  is  assumed  here  that  concentration  polarization  is  negligible  at  electrode 
II.  However,  at  electrode  I,  it  is  significant.  As  a  result,  there  is  a  variation  of 
/itOo  (and  po2 )  through  the  porous  cathode.  The  variation  in  po 2  is  expected 
to  be  approximately  linear.  In  such  a  case,  if  pl0i  p q9  ,  the  net  flux  arriv¬ 
ing  at  the  electrode  I/electrolyte  interface  is  essentially  fixed  (concentration 
polarization  limit),  and  is  proportional  to  (/?q7  —  plQo )  ~  p q2  ,  beyond  some 
threshold  applied  voltage.  In  such  a  case,  the  current  due  to  the  cathodic  re¬ 
action,  namely  5O2  +  2er  — >  O2-,  is  constant,  independent  of  the  applied 
voltage.  However,  the  /Zqo  continues  to  decrease  with  increasing  applied 

voltage.  If  /Zqo  decreases  below  /X^comp,  the  electrolyte  decomposes  result¬ 
ing  in  additional  (excess)  current.  If  the  material  exhibits  nonstoichiometry, 
current  due  to  the  development  of  nonstoichiometry  can  appear  before  de¬ 
composition.  Either  way,  the  measured  current  includes  both  the  cathodic 
current  and  that  due  to  decomposition/nonstoichiometry. 


voltage,  jpqo  <<C  p .  In  such  a  case,  the  cathodic  current  (due 
to  the  cathodic  reaction)  is  fixed.  However,  if  decomposition 
or  nonstoichiometry  develops,  the  measured  current  will  con¬ 
tinue  to  increase  with  the  applied  voltage.  In  such  a  case,  the 
measured  current  density,  /,  versus  the  measured  overpoten¬ 
tial,  r)m,  would  have  little  to  do  with  the  actual  cathodic  re¬ 
action.  Fig.  20  shows  a  schematic  illustration  of  the  cathodic 
current,  the  current  due  to  decomposition/nonstoichiometry, 
and  the  measured  current  as  a  function  of  the  measured  over¬ 
potential.  From  Eq.  (78),  the  minimum  applied  voltage  nec¬ 
essary  to  effect  decomposition  is  given  by 


X  decomp 

II  I  /  ^02  —  ^02 
W  —  m ^  - - - - - 

P  P  \  Ae 


ipi  +  r\  +  r" 
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r o2  ~  Mq2 

\e 


(94) 


Since 


ipi  +  r[  +  r'l 

A 


>  1 


it  is  clear  that  Acp  =  (pu  —  cpl,  the  applied  voltage  necessary 
to  cause  decomposition  given  by  Eq.  (94),  is  lower  than  that 
given  by  Eq.  (92).  That  is,  if  significant  concentration  polar¬ 
ization  exists  at  the  working  electrode,  the  decomposition  (or 
the  development  of  nonstoichiometry)  will  occur  at  a  lower 


Measured  Overpotential 

Fig.  20.  A  schematic  illustration  of  the  measured  current  ( _ )  vs.  the  mea¬ 

sured  overpotential  using  the  three-electrode  system  with  significant  con¬ 
centration  polarization  present.  The  measured  current  consists  of  two  con¬ 
tributions:  (a)  the  cathodic  current  ( _ ),  namely,  ^02  +  2er  — >  O2-,  which 

becomes  concentration  polarization  limited  beyond  some  applied  voltage,  (b) 

The  decomposition  of  nonstoichiometry  current  ( _ ),  which  begins  beyond 

some  applied  voltage.  The  schematic  shows  both  current  and  overpotential 
plotted  on  linear  scales.  If  the  current  is  plotted  on  a  logarithmic  scale,  the 
curves  will  be  convex  up. 

voltage.  If,  however,  r\  r"  and  r\  ^  ipi,  then  Eq.  (94)  re¬ 
duces  to  Eq.  (93).  Several  cases  can  be  examined  correspond¬ 
ing  to  various  relative  values  of  transport  parameters,  and  the 
details  will  albeit  vary.  However,  the  important  conclusion  is 
the  same,  namely  that  in  the  three-electrode  system  under  an 
applied  DC  voltage,  current  due  either  to  the  development  of 
nonstoichiometry  in  the  electrolyte  or  due  to  electrolyte  de¬ 
composition  can  occur,  in  addition  to  the  current  due  to  the  ca¬ 
thodic  reaction.  The  intended  purpose  of  the  three-electrode 
system  is  to  investigate  electrode  kinetics,  i.e.,  to  investigate 
relationship  between  cathodic  overpotential  and  the  current 
due  to  the  cathodic  reaction,  namely  ^02  +  2e/  02~.  The 

net  measured  current,  however,  will  most  likely  include  two 
contributions:  current  due  to  the  cathodic  reaction  and  current 
due  to  the  occurrence  of  decomposition/nonstoichiometry.  It 
is  possible  that  under  certain  conditions  the  current  due  to  de¬ 
composition/nonstoichiometry  may  completely  overwhelm 
the  cathodic  current.  In  such  cases,  the  three-electrode  sys¬ 
tem  will  result  in  significant  overestimation  of  cathodic  elec- 
trocatalytic  activity,  and  may  not  be  a  reliable  technique, 
especially  at  high,  applied  voltages,  and  when  the  mate¬ 
rial  readily  develops  nonstoichiometry — such  as  ceria.  It  can 
be  shown  experimentally  that  current  due  to  nonstoichiom¬ 
etry/electrolyte  decomposition  can  be  several  times  that  due 
to  the  cathodic  reaction  under  relatively  modest  applied  volt¬ 
ages  [22]. 


3.  Discussion 

The  main  objective  of  this  manuscript  was  to  develop 
transport  equations  through  mixed  oxygen  ion  and  electronic 
conducting  membranes,  and  through  predominantly  oxygen 
ion  conducting  membranes,  by  incorporating  transport  across 
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interfaces.  The  analysis  assumes  the  usual  criteria  of  local 
equilibrium  and  electroneutrality.  As  emphasized,  the  as¬ 
sumption  of  local  equilibrium  implies  that  electronic  cur¬ 
rent  cannot  be  assumed  to  be  identically  zero,22  even  in  a 
predominantly  oxygen  ion  conductor,  such  as  zirconia.  This 
has  important  implications  concerning  the  chemical  poten¬ 
tial  of  oxygen  in  the  membrane,  and  thus  the  very  stability 
of  the  membrane.  It  is  well  known  that  many  factors  dic¬ 
tate  transport  through  membranes  and  across  interfaces,  and 
the  assumption  of  transport  parameters  as  constants  (inde¬ 
pendent  of  oxygen  partial  pressure)  is  an  oversimplification. 
At  the  same  time,  however,  such  a  simplification  is  deemed 
necessary  in  order  to  elucidate  the  fundamental  phenomena 
in  a  simple  analytical  framework,  which  is  not  possible  if 
one  incorporates  the  dependence  of  transport  parameters  on 
defect  chemistry  and  on  local  po2 ,23  If  the  exact  dependen¬ 
cies  of  transport  parameters  on  po2  (or  po2 )  are  known  for 
a  given  material  and  interfacial  regions,  which  is  rarely  the 
case,  numerical  solutions  to  transport  equations  can  be  read¬ 
ily  developed.  Such  an  approach  was  not  chosen  here,  as  the 
purpose  was  to  focus  on  fundamental  aspects  of  transport, 
with  emphasis  on  interfaces,  and  not  on  any  specific  material 
with  any  particular  defect  structure.  In  this  framework,  three 
cases  were  examined:  (1)  MIEC  oxygen  separation,  (2)  fuel 
cells,  and  (3)  voltage-driven  oxygen  separation. 

For  the  case  of  MIEC  oxygen  separation,  it  was  shown 
that  the  po2  always  varies  through  the  membrane  monoton- 
ically  between  the  two  end  values  (corresponding  to  the  two 
atmospheres),  regardless  of  the  nature  and  the  details  of  inter¬ 
facial  processes.  This  has  important  implications  concerning 
the  thermodynamic  stability  of  the  membrane.  The  analysis 
shows  that  to  ensure  membrane  stability,  it  is  necessary  that 

^02  >  /U.Q  Comp,  a  condition  which  could  be  achieved  even  if 

^02  <  /XQ2Comp.  Thus,  the  implication  is  that  if  a  membrane 
is  used  for  the  production  of  syngas,  for  example,  it  is  nec¬ 
essary  to  ensure  that  /Xqo  >  /XQ2Comp.  That  is,  the  transport 
properties,  especially  across  interfaces,  must  be  so  adjusted 

to  ensure  that  /Zq2  >  /XQ2Comp.  The  exact  variation  of  \iq2 
through  the  membrane  depends  upon  the  various  transport  pa¬ 
rameters,  including  those  of  interfaces.  This  implies  that  one 
may  need  to  deposit  a  thin  (a  few  nanometers)  layer  of  some 
material  of  appropriate  transport  properties  on  the  permeate 
side  to  ensure  stability  of  an  otherwise  unstable  membrane 
in  the  imposed  atmosphere.  The  same  general  conclusions 
apply  to  fuel  cells.  That  is,  depending  upon  electrode  ef- 


22  Unless  shorted  externally. 

23  There  are  many  published  reports,  which  have  incorporated  the  defect 
chemistry  details  into  the  analysis.  However,  this  necessitated  making  sim¬ 
plifying  assumptions  regarding  the  thermodynamics  of  defects.  Also,  in  so 
doing,  the  important  features  of  electron  transport  across  interfaces  were  gen¬ 
erally  overlooked.  This  not  only  leads  to  incomplete  description  of  transport, 
but  also  does  not  permit  the  estimation  of  /xo,  inside  the  membrane,  unless 
further  simplifying  assumptions  are  made  about  /iq-,  variation  across  inter¬ 
faces.  For  example,  Eq.  (31)  can  explicitly  give  iio2(x)  when  it  is  assumed 
that  /Xq2  =  /Xqo  and  /x'Qi  =  /Zq2  ,  an  assumption  which  is  rarely  justified. 


fects,  an  otherwise  unstable  electrolyte  may  exhibit  stability 
in  fuel.  It  is  known,  for  example,  that  a  fuel  cell  made  with 
phase-stabilized  cubic  or  rhombohedral  Bi203  readily  de¬ 
composes  to  metallic  bismuth  in  hydrogen  as  fuel.  However, 
it  has  been  observed  that  the  same  fuel  cell  remains  stable 
when  operated  in  lightly  humidified  (~3%  H2O)  methane, 

even  when  the  /x o2  in  the  fuel  is  lower  than  /XQ2Comp  [23]. 
This  suggests  that  electrode  characteristics  must  be  such  that 


with  methane  as  a  fuel,  /Xqo  must  be  greater  than  /XQ2comp, 

that  is  /Xqo  >  /Xq2  ,  even  when  /Xqo  <  pQn  .  A  possi¬ 
ble  approach  to  preventing  decomposition  of  bismuth  oxide 
based  electrolyte  is  the  deposition  of  a  more  stable  material 
such  as  YSZ  or  rare  earth  oxide-doped  ceria,  such  as  samaria- 
doped  ceria  (SDC),  on  the  side  exposed  to  fuel.  By  suitably 
tailoring  the  properties  and  thicknesses  of  the  two  regions, 
analysis  based  on  bulk  transport  suggests  that  it  may  be  pos¬ 
sible  to  prevent  electrolyte  decomposition  by  ensuring  that 
the  fio2  in  bismuth  oxide  never  goes  below  the  correspond¬ 
ing  /xJcomp  [14].  However,  the  required  thickness  of  YSZ 
or  SDC  to  prevent  decomposition  is  often  comparable  to  or 
may  even  be  much  greater  than  that  of  bismuth  oxide,  which 
defeats  the  very  purpose  of  using  a  high  conductivity  bismuth 
oxide  based  electrolyte  [14,24].  The  present  work  suggests 
that  a  more  profitable  approach  may  involve  modification  of 
the  anode/electrolyte  interface.  This,  for  example,  could  be 
the  case  III  from  Fig.  7,  where  r"e  ^  re,  r'e.  Thus,  the  present 
work  suggests  a  possible  approach  for  using  BUO3 -based 
electrolytes  by  suitably  modifying  the  anode/electrolyte  in¬ 
terface. 

For  fuel  cells,  equations  for  electrode  overpotentials  in 
terms  of  ionic  and  electronic  current  densities,  and  interfacial 
transport  parameters  were  developed.  The  usual  definition  of 
overpotential  is  as  a  measure  of  the  loss  of  useful  voltage 
in  irreversible  processes.  From  the  standpoint  of  thermody¬ 
namics,  however,  the  overpotential  loss  should  strictly  be  de¬ 
scribed  in  terms  of  the  rate  of  loss  of  potential  work  (which 
can  be  derived  as  decrease  in  free  energy)  into  heat.  The  usual 
approach  is  to  describe  the  overpotential  as  a  measurable  volt¬ 
age  loss,  as  a  function  of  current  (measured  in  the  external 
circuit).  The  inclusion  of  internal  electronic  leakage,  when 
the  overpotential  loss  is  given  as  a  measure  of  the  potential 
work  degraded  as  heat,  leads  to  the  conclusion  that  at  OCV, 
the  overpotential,  defined  as  the  rate  of  loss  of  potential  work 
at  an  electrode  divided  by  the  measured  current,  is  actually 
infinite.  This  is  because  even  when  no  current  is  flowing  in 
the  external  circuit,  part  of  the  potential  work  is  continually 
being  degraded  as  heat  due  to  internal  electronic  leakage. 
As  stated  earlier,  the  electronic  current  can  never  be  zero 
when  potential  gradients  exist,  with  the  result  that  at  OCV, 
the  as-defined  overpotential  is  singular — or  infinite.  Also,  the 
overpotential  as-defined  may  not  be  measurable.  At  the  same 
time,  however,  the  relevant  physical  parameter  is  the  product 
of  the  as-defined  overpotential  and  current,  which  is  always 
finite,  and  is  a  measure  of  the  rate  of  degradation  of  potential 
work  into  heat  (an  irreversible  process).  In  a  typical  practical 
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case  wherein  the  electronic  leakage  is  usually  negligible,  and 
a  finite  current  flows  through  the  external  circuit  (such  that 
\IL\  &  \It\  |4|),  the  overpotential  defined  here  is  essentially 
the  same  as  the  usual  one,  and  measured  experimentally. 

In  both  MIEC  oxygen  separation  membranes  and  fuel 
cells,  the  fio2  and  (p  are  interrelated.  In  voltage-driven  oxy¬ 
gen  separation  systems,  however,  both  /xo2  and  (p  can  be  inde¬ 
pendently  varied,  which  implies  an  additional  (experimental) 
degree  of  freedom.  In  such  a  case,  the  fio2  in  the  membrane 
can  lie  outside  of  the  end  values.  That  is,  /XQ^mbrane,  can  be 
greater  or  lower  than  /Zq2  and/or  /Xq0.  Transport  parame¬ 
ters  of  the  membrane  and  across  interfaces,  and  the  applied 
voltage  determine  the  actual  variation  of  /xo2  in  the  mem¬ 
brane.  Figs.  12-15  show  the  various  possibilities.  An  impor¬ 
tant  point  to  note  is  that  by  virtue  of  the  fact  that  ^j^mbrane 

can  lie  outside  of  the  end  values  implies  that  it  is  possible  for 
^membrane  t0  ke  iower  than  ^^comp #  Under  such  conditions, 

localized  membrane  decomposition  or  nonstoichiometry  can 
occur.  This  has  important  implications  concerning  the  use 
of  the  so-called  three-electrode  system  under  an  applied  DC 
voltage  for  the  measurement  of  electrode  polarization  behav¬ 
ior.  The  analysis  shows  that  in  the  three-electrode  system,  a 
significant  contribution  to  the  net  measured  current  can  oc¬ 
cur  from  membrane  decomposition  (or  the  development  of 
nonstoichiometry),  leading  to  an  inadvertent  overestimation 
of  electrode  electrocatalytic  activity.  The  present  work  thus 
shows  that  caution  should  be  exercised  in  using  the  three- 
electrode  system  under  an  applied  DC  bias  for  the  investi¬ 
gation  of  electrocatalytic  activity  of  an  electrode.  It  is  also 
possible  that  the  /xo2  in  the  membrane  may  be  much  greater 
than  both  /Zq2  and  /Xq?  under  an  applied  DC  voltage.  In  such 
a  case,  cracking  of  the  membrane  can  occur  under  an  applied 
DC  voltage,  as  has  been  previously  shown  theoretically  and 
experimentally  [15,16]. 

All  transport  equations  are  given  under  the  simplifying 
assumption  that  transport  properties  are  independent  of  /jlo2 , 
which  allows  one  to  obtain  simple  analytical  equations.  At 
the  same  time,  the  fact  that  the  chemical  potential  of  oxygen 
is  a  function  of  position,  fio2  (v),  implies  that  small  changes  in 
composition  are  presumed  to  occur  until  steady  state  is  estab¬ 
lished.  The  net  integrated  oxygen  flux  arriving  or  leaving  an 
elemental  region  dx  is  directly  related  to  the  interrelationship 
between  composition  (stoichiometry)  and  fio2 .  If  for  a  given 
change  in  /xo2,  namely  A/xo2,  the  change  in  composition  is 
small  (narrow  stoichiometry  range),  this  means  the  capacity 
of  the  internal  EMF  will  also  be  small.  If  chemical  capacitors 
are  used  to  describe  the  equivalent  circuit,  then  it  will  im¬ 
ply  that  the  magnitude  of  the  chemical  capacitance  is  small. 
This  may  likely  be  the  case  with  zirconia.  Alternatively,  if 
for  a  given  change  in  /xo2 ,  namely  A/xo2 ,  the  change  in  com¬ 
position  is  large  (wide  stoichiometry  range),  this  means  the 
capacity  of  the  internal  EMF  will  also  be  large.  If  chemical 
capacitors  are  used,  then  it  will  imply  that  the  magnitude  of 
the  chemical  capacitance  is  large.  This  may  likely  be  the  case 
with  ceria.  Regardless  of  whether  internal  EMF  or  chemical 
capacitors  are  used,  it  is  clear  that  a  complete  description  of 


transport  will  require  relationship  between  composition  (sto¬ 
ichiometry)  and  fio2 ,  which  will  most  likely  have  to  be  de¬ 
termined  experimentally.  Unfortunately,  this  is  rarely  known 
for  any  material  with  any  degree  of  certainty. 

4.  Summary 

In  the  present  manuscript  equations  were  developed  to  de¬ 
scribe  transport  of  oxygen  ions  and  electrons  through  mem¬ 
branes  by  explicitly  taking  into  account  interface  properties. 
The  fundamental  assumptions  made  are  the  same  as  in  virtu¬ 
ally  all  of  the  earlier  studies,  namely,  the  existence  of  elec¬ 
troneutrality  and  local  equilibrium.  The  assumption  of  local 
equilibrium  is  shown  to  imply  that  electronic  current  can¬ 
not  be  entirely  neglected,  however  small,  to  ensure  that  lo¬ 
cal  chemical  potential  of  oxygen,  fio2 ,  is  uniquely  defined. 
That  is,  even  in  a  predominantly  oxygen  ion  conductor,  elec¬ 
tronic  transport  must  be  taken  into  account.  Transport  of 
oxygen  ions  and  electrons  through  the  membrane  as  well 
as  across  electrode/membrane  interfaces  was  taken  into  ac¬ 
count  in  the  present  work.  These  equations  were  then  applied 
to  fuel  cells,  MIEC  oxygen  separation  membranes  under  an 
applied  pressure  gradient,  voltage-driven  oxygen  separation, 
and  the  three-electrode  system  which  is  often  used  for  the 
study  of  electrode  kinetics.  It  was  shown  that  in  fuel  cells  or 
MIEC  oxygen  separation  membranes,  as  long  as  the  chemi¬ 
cal  potential  of  oxygen  on  the  reducing  side  (anode  in  the 
case  of  fuel  cells;  permeate  side  for  MIEC  oxygen  sepa¬ 
ration  or  fuel  side  for  the  synthesis  of  syngas  using  MIEC 
membranes)  is  greater  than  the  chemical  potential  of  oxygen 
corresponding  to  membrane  thermodynamic  stability,  mem¬ 
brane  decomposition  should  not  occur.  Membrane  may  con¬ 
tinue  to  be  stable  even  if  the  chemical  potential  of  oxygen 
in  the  anode  gas  (fuel  cells)  or  permeate  side  (MIEC  mem¬ 
branes)  is  lower  than  the  decomposition  chemical  potential, 
provided  the  interface  transport  properties  are  suitably  tai¬ 
lored  such  that  the  chemical  potential  of  oxygen  just  inside 
the  interface  is  greater  than  the  decomposition  chemical  po¬ 
tential,  that  is  /Xqo  >  /XQ2Comp.  However,  in  the  case  of  oxy¬ 
gen  separation  under  the  action  of  an  applied  voltage,  de¬ 
pending  upon  the  magnitude  of  the  applied  voltage  and  the 
relative  transport  parameters  of  the  membrane  and  interfa¬ 
cial  regions,  the  chemical  potential  of  oxygen  in  the  mem¬ 
brane  inside  the  interfaces  can  go  below  the  decomposition 
chemical  potential.  In  such  a  case,  membrane  decomposi¬ 
tion  can  occur,  or  a  nonstoichiometry  may  develop.  This  has 
significant  implications  concerning  the  use  of  the  so-called 
three-electrode  system  for  the  study  of  electrode  kinetics,  for 
example  for  the  study  of  a  prospective  cathode  for  SOFC. 
The  significance  is  that  measurements  made  using  the  three- 
electrode  system  under  high  applied  voltages  may  yield  er¬ 
roneous  conclusions  concerning  cathode  activity,  since  the 
net  current  generated  will  consist  of  two  terms:  (a)  that  due 
to  the  actual  cathodic  reaction,  namely  ^(4  +  2er  — >  O2-, 
and  (b)  that  due  to  the  occurrence  of  decomposition  and/or 
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the  development  of  nonstoichiometry  (oxygen  deficiency).  In 
a  typical  experiment,  the  measured  current  is  invariably  at¬ 
tributed  to  the  cathodic  reaction,  thereby  overestimating  the 
cathodic  activity  of  the  prospective  cathode,  since  the  exper¬ 
imenter  has  no  way  of  separating  the  two.  In  many  cases, 
the  current  due  to  nonstoichiometry/decomposition  may  be 
much  larger  than  the  cathodic  current,  leading  to  large  er¬ 
rors  in  the  estimation  of  cathodic  activity.  The  analysis  also 
shows  that  the  decomposition  (or  the  development  of  oxy¬ 
gen  deficiency)  may  occur  just  under  the  positive  electrode 
or  just  under  the  negative  electrode,  depending  upon  the  rel¬ 
ative  magnitudes  of  ionic  and  electronic  transport  parame¬ 
ters  of  the  bulk  and  the  interfaces.  In  practice,  however,  the 
occurrence  of  decomposition  (or  the  development  of  oxy¬ 
gen  deficiency)  under  the  positive  electrode  is  deemed  un¬ 
likely.  This  is  because  in  such  a  case,  the  electrons  required 
for  the  neutralization  of  the  cation  (Zr4+  to  Zr)  transport 
through  the  membrane,  which  may  be  kinetically  hindered 
if  the  electronic  conductivity  is  low.  Decomposition  (or  the 
development  of  oxygen  deficiency)  under  the  positive  elec¬ 
trode,  however,  is  likely  in  the  case  of  ceria-based  mem¬ 
branes  under  appropriate  conditions  since  electronic  con¬ 
ductivity  of  ceria  is  typically  much  higher,  especially  at  low 

PO2  • 

All  transport  equations  are  given  in  a  general  form,  with¬ 
out  invoking  any  particular  defect  structure.  The  rationale  is 
that  fundamental  stability  issues  are  not  dependent  on  the 
types  and  nature  of  defects.  All  transport  equations  are  given 
assuming  transport  parameters  to  be  independent  of  fi o2 .  Be¬ 
cause  of  this  assumption,  the  spatial  variations  of  \io2  >  and 
/Zq2—  through  the  membrane  are  shown  to  be  linear.  This  as¬ 
sumption  also  allows  the  use  of  simple  equivalent  circuits. 
Bulk  regions  were  modeled  using  a  combination  of  ionic  and 
electronic  area  specific  resistances  for  the  bulk  and  an  inter¬ 
nal  EMF  source.  Interfacial  regions  were  also  modeled  using 
internal  EMF  sources,  instead  of  capacitors,  in  combination 
with  ionic  and  electronic  area  specific  resistances.  This  is 
entirely  consistent  with  the  definition  of  Nernst  voltage  in 
terms  of  the  net  A/xo2  across  an  interface.  As  such,  its  mod¬ 
eling  as  an  EMF  source  is  equally  appropriate,  and  in  fact,  is 
more  consistent  with  the  general  thermodynamic  framework. 
Usually  the  transport  parameters  are  /zo2  -dependent.  In  such 
cases,  the  variation  of  /xo2 ,  <p,  and  jl02-  will  not  be  spatially 
linear.  If  the  dependencies  for  a  given  material  and  interfa¬ 
cial  regions  are  precisely  known,  which  is  rarely  the  case,  the 


spatial  variations  of  /xo2,  <p,  and  jl02-  can  be  easily  evalu¬ 
ated  numerically — if  not  analytically.  Such  an  approach  was 
not  selected  here.  The  simplification  made  may  only  alter  the 
details,  but  not  the  broader,  general  conclusions. 
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